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ABSTRACT 

t~^ ' The giant merger remnant galaxy NGC 1316 (Fornax A) is an ideal probe for studying the 

CsJ \ long-term effects of a past major merger on star cluster systems, given its spectroscopically 

derived merger age of ^ 3 Gyr which we reported in a recent paper (Goudfrooij et al. 2001). 
Here we present new ground-based, large-area optical and near-IR imaging of star clusters in 
NGC 1316. The ground-based photometry is complemented with deep Hubble Space Tele- 
scope WFPC2 imaging, constituting an excellent combination for studying globular cluster 
systems. We find that the optical-near-IR colours and luminosities of the brightest ^ 10 clus- 
ters in NGC 1316 are consistent with those of intermediate-age (2-3 Gyr) populations. In 
f^ ■ particular, the near-IR data preclude ages ^ 1.5 Gyr and ^ 4 Gyr for those clusters. Un- 

like 'normal' giant ellipticals, the B — I colour distribution of clusters in NGC 1316 is not 
clearly bimodal. However, the luminosity functions (LPs) of the blue and red parts of the 
cluster colour distribution are different: The red cluster LF is well represented by a power 
law, (j){L) dL (X 2^-i-2±o.3 ^^^^ extending to about 1.5 mag brighter (in B) than those of typ- 
ical giant ellipticals. In contrast, the shape of the blue cluster LF is consistent with that of 
O ' 'normal' spiral and elliptical galaxies. We conclude that the star cluster system of NGC 1316 

is a combination of a population of age ^ 3 Gyr having roughly solar metallicity and a popu- 
^ ■ lation of old, metal-poor clusters which probably belonged to the pre-merger galaxies. After 

^ ' the 3 Gyr old, metal-rich clusters fade to an age of 10 Gyr, they will form a red 'peak' in 

a bimodal cluster colour distribution. This 'red peak' will have a colour consistent with that 
Lj found in 'normal, old' giant ellipticals of the same galaxy luminosity (taking age dimming 

rS I into account). The surface density profile of clusters in the innermost regions is lower than 

jrt ■ that of the integrated light of the galaxy, presumably due to the collective effect of extended 

star formation in the inner regions of NGC 1316 and tidal shocking of the inner clusters. Out- 
side the core, the surface density profile of clusters is consistent with that of the underlying 
starlight, suggesting that the cluster system originally experienced the same violent relaxation 
as did the main body of the merger remnant. The specific cluster frequency is presently Sn 
= 1 .7 ± 0.4 down to the 50 per cent completeness limit of the WFPC2 photometry, and will 
increase to Sn ^ 2.0 as the merger-induced stellar (and star cluster) population fades to an 
age of '^ 10 Gyr (barring further merger events), consistent with specific frequencies of typ- 
ical giant ellipticals in the field and in poor groups. These features of the star cluster system 
of NGC 1316 are fully consistent with scenarios for forming 'normal' giant elliptical galaxies 
through gas-rich mergers at look-back times <; 10 Gyr 

Keywords: galaxies: individual: NGC 1316-galaxies: elliptical -galaxies: radio -galaxies: 
interactions - globular clusters: general 



1 INTRODUCTION 

* E-mail (internet): goudfroo@stsci.edu 

t Affiliated to the Astrophysics Division, Space Science Department, Eu- Globular clusters (hereafter GCs) are among the few observable 

ropean Space Agency fossil records of the formation and evolution of galaxies. Their na- 
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ture as simple stellar population (hereafter SSP) significantly sim- 
plifies the determination of their ages and metallicities relative to 
that of the stellar populations that constitute the integrated light of 
their parent galaxies. In our Galaxy, they have yielded crucial in- 
formation on early chemical enrichment, the formation time scale 
of the galactic halo, the presence of a 'thick disk', and an estimate 
for the age of the universe. While the relative faintness of extra- 
galactic GCs has long prevented one from obtaining high-quality 
photometry and (especially) spectroscopy, the Hubble Space Tele- 
scope (HST) and large ground-based telescopes now allow this goal 
to be achieved. 

One important, well-known aspect of GC systems among 
galaxies is that the number of clusters per unit galaxy luminosity 
(named the specific frequency Sn) increases systematically from 
late-type to early-type galaxies, being ~ 2-3 times higher in ellip- 
tical (E) than in spiral galaxies of type Sb and later (Harris & van 
den Bergh 1981; Harris 1991). This fact has been used as an argu- 
ment against the scenario for forming ellipticals through mergers: 
van den Bergh (1995) argued that if ellipticals form through gas- 
rich mergers during which star (and GC) formation occurs with 
a normal (e.g., Salpeter (1955)) initial mass function (IMF), the 
^jv of the resulting merger remnant galaxies should not be signif- 
icantly different from the progenitor galaxies since both GCs and 
stars would have been formed during the merger. On the other hand, 
recent observations with the HST have led to a wealth of discover- 
ies of young GCs in merging and starburst galaxies (e.g., Holtz- 
man et al. 1992; Whitmore et al. 1993, 1999; Meurer et al. 1995; 
Schweizer et al. 1996; Miller et al. 1997). Indeed, the formation 
of young stars in high-density regions of starbursts seems to occur 
preferentially within star clusters rather than unbound associations 
(Meurer et al. 1995; Schweizer et al. 1996), most probably from 
giant molecular clouds whose collapse is being triggered by a 100- 
to 1000-fold increase in gas pressure due to supernova and shock 
heating during starbursts (e.g.. Jog & Solomon 1992; Elmegreen & 
Efremov 1997). As gas-rich galactic mergers produce the most en- 
ergetic starbursts known (e.g., Sanders & Mirabel 1996), the higher 
specific frequency of GCs in ellipticals relative to that in spirals 
may be (partly) accounted for by secondary populations of GCs 
created during gas-rich mergers. 

In this respect, a particularly interesting feature of the GC sys- 
tems of many giant ellipticals is the presence of bimodal colour 
distributions, providing clear evidence for the occurrence of a 'sec- 
ond event' in the formation of these systems (see, e.g., the review 
by Ashman & Zepf 1998 and references therein). While such a bi- 
modal colour distribution was actually predicted from the merger 
scenario of galaxy formation (e.g.. Ashman & Zepf 1992), opinions 
about the general nature of the 'second event' differ among authors, 
as detailed upon below. One important constraint to any scenario to 
explain the bimodality is set by the colour of the 'red' peak in the 
colour distribution of GCs among giant ellipticals, which (if inter- 
preted in terms of metallicity) correlates well with the metallicity 
and luminosity of the parent galaxy. In contrast, the mean metal- 
licity of the 'blue' peak is more or less constant among galaxies 
as reported by several groups (Forbes, Brodie & Grillmair 1997; 
Forbes & Forte 2001; Kundu & Whitmore 2001; but see Larsen et 
al. 2001 for a somewhat different view). There are three compet- 
ing scenarios that attempt to explain the bimodality of GC colours 
in giant ellipticals. The 'merger model' (Schweizer 1987; Ashman 
& Zepf 1992) and the 'multi-phase collapse model' (Forbes et al. 
1997) both suggest that it is the consequence of GCs forming dur- 
ing two distinct epochs in the star formation history of these galax- 
ies, whereas the 'accretion scenario' (Forte, Martinez & Muzzio 



1982; Cote, Marzke & West 1998) proposes that every galaxy is 
bom with a GC system that has a median colour (metallicity) ac- 
cording to the colour-magnitude relation among galaxies (e.g., Ter- 
levich et al. 1999), and that the bimodality is due to accretion of 
small galaxies (with their associated metal-poor GCs) by a larger 
galaxy (with its pre-existing metal-rich GCs). The first two models 
differ in the formation mechanism of the metal-rich clusters: the 
'merger model' suggests that they are formed during a major, gas- 
rich galaxy merger, while the 'multi-phase collapse model' pro- 
poses that they are formed during a secondary collapse phase of 
the galaxy, which does not involve a merger. 

Given these different points of view, it is important to obtain 
any additional evidence for (or against) the scenario in which the 
second-generation GCs are formed during gas-rich galaxy merg- 
ers. In this context it is of great interest to find out what the ages 
and metallicities of the luminous GCs in merger remnants actu- 
ally are, and whether or not the photometric properties of the GC 
populations found in merger remnants will evolve into those typ- 
ically found in 'normal' giant ellipticals. The currently most de- 
pendable age determinations stem from Hydrogen Balmer lines ob- 
served in spectra of one bright GC in the peculiar elliptical galaxy 
NGC 1275 (Zepf e t al. 1995a; Brodie et al. 19 98) and of two such 
GCs in NGC 7252 ( [Schweizer & Seitzer 199^ ). The luminous GCs 
in these merger remnants were found to be ~ 0.5 Gyr old (cf. also 
Hibbard & Mihos 1995; Schweizer 1998). Unfortunately however, 
the metallicities of these clusters were hard to measure accurately, 
since metallic lines are intrinsically weak at this age (the spectrum 
being dominated by A-type stars). Thus, it still remained to be seen 
whether or not those luminous clusters could evolve into the 'red' 
clusters seen in 'normal' giant ellipticals. Another important test 
of the 'merger scenario' is to find out whether or not these clus- 
ters survive the merger era, i.e., whether or not they can be found 
in merger remnants at different times after a (major) merger. This 
issue has been brought up for the case of the ~ 0.5 Gyr old clusters 
in NGC 1275: Brodie et al. (1998) found H7 and US equivalent 
widths that are somewhat larger (and B~R colours that are bluer) 
than those of the Bruzual & Chariot models that they used, both for 
Salpeter or Scalo (1986) IMFs. They showed that the large Balmer 
line equivalent widths could be brought into agreement with those 
models by assuming a truncated IMF involving a mass range of 2-3 
M0. For such a high low-mass cutoff, they argue that these clusters 
should fade away in only ~ 1 10^ yr. On the other hand, the ^ 0.5 
Gyr old clusters in NGC 7252 show Balmer line strengths that are 
well reproduced by the Bruzual & Chariot (1996) models (hereafter 
the BC96 models) that use the Salpeter or Scalo IMF (Schweizer 
& Seitzer 1998). Also, Gallagher & Smith (1999) found no need 
to invoke a truncated IMF to explain the Balmer line strengths in 
young clusters in the starburst galaxy M 82. But then again, a de- 
tailed study of the very luminous star cluster M 82-F (Smith & 
Gallagher 2001) revealed that a lower mass cutoff of 2-3 Mq is 
required to match their observations for a Salpeter IMF. Obviously, 
the IMF among luminous young star clusters may not be universal. 

We (Goudfrooij et al. 2001, hereafter Paper I) recently ob- 
tained multi-slit spectra of GC candidates in NGC 1316 (Fornax 
A), an early-type merger remnant galaxy (see, e.g., Schweizer 
1980). We discovered the presence of ~ 10 GCs associated with 
NGC 1316 that have luminosities up to an order of magnitude 
higher than that of uj Cen, the most massive cluster in our Galaxy. 
Our measurements of Ha and the Ca II triplet in the spectra of the 
brightest GCs showed them to have solar metallicity (to within ± 
0.15 dex) and to have an age of 3.0 ± 0.5 Gyr. This reinforces 
the view that luminous GCs such as those found before in younger 
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merger remnants do not necessarily have high low-mass cutoffs to 
their IMFs, and it means that they can actually survive disruption 
processes taking place during a galactic merger 

In the present paper we study the photometric properties of 
NGC 1316's GC system. We use a combination of HST WFPC2 
photometry and ground-based, large-field optical and near-IR im- 
agery of NGC 1316. The uniquely high spatial resolution of the 
WFPC2 data provides accurate magnitudes and colours of the GCs, 
with negligible contamination by foreground stars or background 
galaxies. The large-field optical data are used to determine accu- 
rate values for the cluster specific frequency, while the near-IR 
data allow us to significantly improve age dating and to search 
for intermediate-age, Magellanic-type GCs wh ich feature very red 
optical-IR colours, such as found in NGC 5 128 ( Minnitietal. 199^ ) 
and NGC 7252 (Maraston et al. 2001). 

This paper is built up as follows. Section eI summarizes main 
properties of NGC 1316. The different observations and data re- 
duction methods are described in section pi The various results on 
the photometric properties of the GC system are presented in Sec- 
tion W. Section pi summarizes the main conclusions of this study. 
Finally, Section^ briefly discusses important outstanding issues. 



2 NGC 1316 

The giant early-type galaxy NGC 1316 (= Fornax A = PKS 
0320 - 37 = ESQ 357 - G022 = ARP 154) is one of the bright- 
est (Bt = 9.42 mag, de Vaucouleurs et al. 1991, hereafter RC3) 
and closest radio galaxies in the sky. It is located in the outskirts of 
the Fornax cluster, at a projected distance of 3? 7 from NGC 1399, 
the central giant elliptical galaxy. Several features of NGC 1316 
establish firmly that it is a merger remnant. The outer envelope in- 
cludes several non-concentric arcs, tails and loops that are most 
likely remnants of tidal perturbations, while the inner part of the 
spheroid is characterized by non-concentric ripples of 0.1-0.2 mag 
amplitude above a best-fitting r^'* law, a surprisingly high cen- 
tral surface brightness and small core radius (and effective radius) 
for its galaxy luminosity (Schweizer 1980, 1981; Caon, Capaccioli 
& D'Onofrio 1994). These characteristics, together with a velocity 
dispersion which is significantly lower than that of other elliptical 
galaxies of similar luminosity, cause NGC 1316 to lie far off the 
fundamental plane of early-type galaxies or the Faber- Jackson re- 
lation (e.g., D'Onofrio et al. 1997). Figure hi shows a B-band image 
of the inner 7!5 x 7^5 of NGC 1316 in which some of the features 
mentioned above can be found. All of these features are consistent 
with NGC 1316 having undergone a galactic merger after which 
dynamical relaxation has not yet had time to complete fully. 

As to the adopted distance of NGC 1316, we take advantage 
of the fact that two well-observed type la supernovae (SNe la) have 
occurred in NGC 1316 (SN1980N and SN1981D). Using the pre- 
cise distance indicator for SNe la that utilizes the tight relation be- 
tween their light curve shape, luminosity, and colour (Riess et al. 
1998 and private communication), we arrive at a distance of 22.9 
(± 0.5) Mpc for NGC 1316, equivalent to (m- A/)o = (31.80 ± 
0.05). At this distance, 1" corresponds to 111 pc. Note that this 
sets NGC 1316 slightly behind the core of the Fornax cluster (for 
which (m-A/)o = 31.54 ± 0.14; Ferrarese et al. 2000). In fact, 
several pieces of evidence suggest that NGC 1316 and a group of 
smaller galaxies around it form a dynamical entity that is separate 
from the main Fornax cluster: (i) its location in the periphery of 
the Fornax cluster, (ii) its luminosity that is 1 . 1 mag brighter than 
that of NGC 1399, the brightest giant elliptical in the Fornax clus- 



ter core, and (Hi) the clustering of the radial velocities of galaxies 
near NGC 1316, with a median radial velocity that is ^ 300 km 
s~^ higher than that of the central regions of the Fornax cluster 
(e.g., Ferguson 1989; RC3). The dynamical study of the Fornax 
cluster by Drinkwater, Gregg & Colless (2001) indeed shows that 
the Fornax system has two components: the main Fornax Cluster 
centered on NGC 1399 and a small subcluster centered on NGC 
1316. Drinkwater et al. show that this partition is preferred over a 
single cluster at the 99% confidence level. We therefore assume in 
the following that the envirormient of NGC 1316 is a poor group 
rather than a rich cluster. Other global properties of NGC 1316 are 
listed in Table 1 of Paper I. 



3 OBSERVATIONS, DATA REDUCTION AND CLUSTER 
PHOTOMETRY 

3.1 HST WFPC2 imaging 

3.1.1 Observations and data reduction 

We retrieved archival images of NGC 1316 taken with the Wide 
Field and Planetary Camera 2 (WFPC2) aboard HST in order to 
study the star clusters in the iimer regions of this galaxy. The data 
consist of multiple images through the F450W and F814W fil- 
ters. The total exposure times were 5000s in F450W and 1860s 
in F814W. A subset of the images were spatially offset by 0'.'5 
from the others, which corresponds to an approximately integer 
pixel shift in both PC and WF CCDs. After standard pipeline pro- 
cessing and alignment of the images, we combined the images us- 
ing the STSDAS task CRREJ, which effectively removed both cos- 
mic rays and hot pixels. We also trimmed each image to exclude 
the obscured regions near the pyramid edges of WFPC2, yielding 
751 X 751 usable pixels per CCD. Results from the authors of the 
HS T prog ram (proposa l ID 5 990) have been publi shed in Forbes et 
al. ( |l998| ), Elson et al. ( |l998| ) and Grillmair et al. ( |l999| ). 

One of the main sources of uncertainty in WFPC2 photometry 
is due to charge-transfer efficiency (CTE) problems of the CCDs 
for which correction recipes are available (Holtzman et al. 1995; 
Whitmore, Heyer & Casertano 1999a). In the present case however, 
compact-source photometry on a high and strongly varying back- 
ground from NGC 1316, we ignored any such correction. The low- 
est background surface brightness was encountered in the F450W 
WF images, but even there all photometered objects lie in regions 
where the background was in excess of 10 counts per pixel. The 
lack of CTE correction is therefore expected to affect the magni- 
tudes by less than 0.02 mag and the colours by less than 0.01 mag, 
negligible for the purposes of this paper 

Prior to performing source photometry, the strongly varying 
galaxy background was fitted. The main reason for this is to min- 
imize errors in the photometry due to any particular choice of ob- 
ject aperture and sky annulus (see below). For the WF images, the 
smooth background gradient was fitted by a bi-cubic spline fit and 
subtracted. Care was taken not to include any point-like source or 
dust filament in the fit. For the PC image, any fit to the galaxy 
background is hampered by the presence of very prominent dust 
features reaching into the nucleus (cf. Fig. 111). After extensive ex- 
perimentation, we found the following procedure to give the best 
results. We first built an isophotal model composed of pure ellipses 
(separately for each filter) using the IRAF/STSDAS task ELLIPSE in 
the ISOPHOTE package. Care was taken to mask out all pixels cov- 
ered by the dust features prior to the fitting procedure, as well as 
bright sources, using judgement by eye. Subtraction of the isopho- 
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Figure 1. Grey-scale plot of _B-band CCD image of the central 7^5 X 7.'5 of NGC 1316, observed at the NTT (see Section |3^ 
by the bar at the lower left of the image. North is up and East is to the left. 



. The plate scale is indicated 



tal model revealed a 'residual' image on which compact objects as 
well as the dust features stand out. The galaxy background in the 
parts of the PC image that were not covered by the elliptical model 
fit was fit by applying a median filter with a 1 1 x 1 1 pixel kernel 
to the PC image. After setting the intensity of the median-filtered 
image to zero in the locations where the isophotal model was non- 
zero, the final 'model image' for the PC was created by summing 
the isophotal model and the median-filtered image together. 



3.1.2 Cluster candidate selection and photometry 

The selection and photometry of star cluster c andidates wa s car- 
ried out using the IRAF version of DAOPHOT-II ( [Stetson 1987 ). The 
objects were selected by applying the DAOFIND task to an image 
prepared by dividing the F450W image by the square root of the 
model image (thus having uniform shot noise characteristics over 
the whole image). We adopted fairly tight shape constraints [—0.6 
< roundness < 0.6; 0.2 < sharpness < 0.9] in order to exclude ex- 
tended background galaxies and faint objects distorted by noise or 
any residual bad pixels. The detection threshold was set at 4 sigma 
above the residual background. Although DAOFIND returned with 
apparent point-like detections located within -or on the edge of- 
the dust features, we decided to exclude those from further analysis 
due to the difficulty in judging the location of those sources relative 
to the dust features along the line of sight, which is bound to yield 
spurious results. 



In order to convert the STMAG magnitudes F450W and 
F814W into ground-based Johnson B and Cousins / magnitudes, 
we used the SYNPHOT package within STSDAS. Synthetic spec- 
tra of ste llar ty pes ranging from BOV to K3III from the Bruzual & 
Chariot (1993 1 library were convolved with Johnson B and Cousins 
/ filters to yield magnitudes relative to Vega, and with F450W and 
F814W filters to yield magnitudes in the STMAG system. Polyno- 
mial fits between the two different magnitude systems (to second 
order in (F450W — F814W)) resulted in the following conversions: 



B - F450W 



/-F814W 



(0.643 ± 0.002) 

+ (0.138 ± 0.005) (F450W 
+ (0.024 ± 0.002) (F450W 

(-1.267 ±0.001) 

+ (0.016 ± 0.003) (F450W 
+ (0.011 ± 0.002) (F450W 



F814W) 
F814W)^ 

F814W) 
F814W)^ 



Aperture photometry was carried out through small apertures 
in order to reduce the noise from the residual background. The 
background level was determined in an annulus surrounding each 
object aperture. After extensive experimentation with different ob- 
ject aperture sizes, we adopted apertures of 2 pixel radius (r) for 
both the PC and WF CCDs. This choice was made after consider- 
ing (i) that larger-aperture radii yielded significantly larger errors 
for derived B — I colours, and (ii) that smaller aperture radii for 
the WF CCDs (i.e., 1 pixel) would be expected to cause system- 
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Figure 2. Difference between B — I colour indices of point-like sources 
measured on the unbinned PC CCD ('PC on the plot) and those measured 
on the PC in 2 X 2 binned form ('BPC on the plot). Solid circles are mea- 
surements on the binned PC using an aperture radius of 2 (binned) pixels, 
and open cicles are measurements on the binned PC using an aperture ra- 
dius of 1 pixel. The measurements on the unbinned PC were made with an 
aperture radius of 2 pixels. Note the relative insensitivity of the photometric 
accuracy to the choice of the aperture radius. 



atic biases in the final photometry such as centering errors and any 
variation of aperture corrections with position on the CCD. On the 
other hand however, small aperture radii would yield superior pho- 
tometric accuracy in the presence of significant variations in the 
background. This is why we tested whether a smaller aperture ra- 
dius ( 1 WF pixel) for the WF CCDs would improve the photomet- 
ric accuracy, as follows. We compared the colours of the point-like 
sources on the PC with that of the same sources after having binned 
the PC CCD by 2 X 2 pixels, to simulate the pixel size of the WF 
CCDs. We performed this comparison with (a) r = 2 pix for the 
unbinned PC and r = 1 pix for the binned PC, and (b) r = 2 pix 
for the unbinned PC and r = 2 pix for the binned PC. A perture cor- 
rections were performed according to Holtzman et al. (1995). We 
depict the results of this comparison in Fig. 



It can be seen that 
the accuracy of the photometry is not very sensitive to the choice of 
aperture radius. The standard deviation of A(B— /), the difference 
of the colours from the binned and unbinned PC, was 0.14 for case 
(a) described above, and 0.11 for case (b). This encouraging result 
is mainly due to the successful removal of the galaxy background 
as described in the previous subsection. 

Galactic foreground ext inction towards NGC 1316 is Ab ~ 
0.0 ( ^urstein & Heiles 1984| )- 

A B vs. B—I colour-magnitude diagram (hereafter CMD) for 
the 492 detected objects is shown in Fig. 0, in which the objects in 
the PC frame are assigned a symbol different from those in the WF 
frames. The direction of the reddening vector is also shown. Note 
that the objects in the PC frame are typically somewhat redder than 
those in the WF frames, as could be expected due to the extensive 
dust filaments in the PC frame. From this list of detections, a strict 
cluster candidate list was made by applying a colour selection cut 
of 1.0 < B — I < 2.5, which is the full range expected for clus- 
ters older than 1 Gyr and -2.5 < [Z/Zq] < 0.5 (BC96; Maraston 
1998; see also Section 4.1 below). The quoted lower limit to the 
age of the clusters can be considered realistic since we found previ- 
ously from spectroscopy that even the brightest GCs in NGC 1316 
have an age of ~ 3 Gyr (Paper I). For objects in the PC frame, the 
upper limit to B — / was chosen to be 3.0, allowing for Ab Ss 1 



m 23 




N"GC 1316 
WTPCS data 



B-I 

Figure 3. B vs. B — I colour-magnitude diagram for NGC 1316 GC can- 
didates from the WFPC2 data. Only objects with B — I eiTors less than 0.5 
mag are plotted. Objects in the PC frame are shown as open circles, while 
objects in the WF frames are shown as solid circles. The dotted lines de- 
hneate the colour and magnitude criteiia for cluster candidates, while the 
hatched region represents the apparent magnitude and colour range for GCs 
in the halo of our Galaxy, placed at the distance of NGC 1316 (see text for 
details). The solid line shows the reddening vector for Ab = 1 mag of 
extinction. 



mag of extinction. The CMD shown in Fig. H is further discussed in 



Section 4. 1 



Photometry and astrometry for the 50 brightest cluster 

in the Ap- 



candidates in the WFPC2 frames is shown in Table Al 



pendix. A similar table for the full list of cluster candidates in the 
WFPC2 frames can be obtained in electronic form from CDS. 

Completeness tests were done by adding artificial GCs to the 
images and re-applying the aperture photometry programs in a way 
identical to those used for the original images. Artificial GCs were 
generated from a composite Point Spread Function (PSF) derived 
from real GC candidates in each frame, and added in batches of 
100 at different magnitude intervals, at random positions within the 
frames, with colours equal to the median colour of GC candidates. 
The resulting completeness function for the WF frames is shown in 
Fig. W. The 50 per cent completeness level is at B ~ 25.6 for the 
WF frames, and B ~ 24.5 for the PC frame. 



3.2 Ground-based optical imaging 

3.2.1 Obser-vations and data reduction 

Ground-based CCD images of NGC 1316 were retrieved from 
the ESO archive database. The data were originally obtained in 
November 1992 at the ESO 3.5-m New Technology Telescope 
(NTT) with the red arm of the ESO Multi-Mode Instrument 
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Figure 4. S and /-band completeness function for the WFPC2 photom- 
etry of NGC 1316 cluster candidates on the WF frames. The 50 per cent 
completeness level is at i? ^^ 25.6 and / ~ 25.2 mag. 



(EMMI). The detector was a thick, front-illuminated CCD of type 
Thomson 31156 Grade A, having 1024x 1024 sensitive pixels. The 
pixel size was 19/im x 19/im, yielding a scale of 0.44 arcsec per 
pixel and a total field of view of 7'5 x 7'5 (see Fig. |i|). NGC 1316 
was observed through B, V and Ic (Cousins) filters, with total ex- 
posure times of 4200 s, 2100 s and 900 s, respectively, split up into 
several exposures. 

After retrieving several bias, dark, and twilight flat field 
frames taken during the observing runs in which the images of 
NGC 1316 were taken, basic data reduction was carried out in the 
standard way using the IRAF CCDRED package. Cosmic radiation 
events in the individual images were removed by an appropriate 
averaging program. First, the images taken through the same filter 
were aligned to a common coordinate system, using the centroids 
of stars in the field of view. This alignment procedure was accurate 
to within 0.03 pixel. The images were subsequently averaged to- 
gether by comparing all individual pixel values (per unit exposure 
time) with the median value over all frames (taken using the same 
filter). Individual pixels are rejected during the combining process 
if their value exceeds the range expected from the (sky -I- read-out) 
noise. The effective seeing values of the combined images are l'.'28 
in B, 1'.'23 in V, and 1" 17 in Ic- Fig. Ill shows the final combined 
B-band image. 

Photometric calibrations were derived using images of the 
Rubin 149 standard star field for which the photometry in the 
Johnson-Kron-Cousins system is given in Landolt (1992). The re- 
sulting colour equations were: 



B-b = 


= (22.977 ± 0.002) + (0.032 ± 0.002) (6- 


-V) 


V -V = 


= (24.466 ± 0.005) + (0.027 ± 0.005) (6- 


-v) 


I-i = 


= (24.060 ± 0.002) - (0.020 ± 0.003) («- 


-i) 



where upper case BVI denote calibrated Johnson-Kron-Cousins 
magnitudes, and lower case bvi denote instrumental magnitudes 
(normalized to an exposure time of 1 s). We used the standard atmo- 
spheric extinction coefficients for the ESO site (ESO User Manual, 
1993). This calibration was double checked with WFPC2 photom- 
etry of stars in the area covered by the WFPC2 images and with 
published aperture photometry of Wegner (1979). The differences 
were minor (less than 0.02 mag in any band). 



3.2.2 Cluster candidate selection and photometry 

As described above for the case of the HST WFPC2 images, we 
attempted to remove the galaxy background from the images. First 
we used the STSDAS task ELLIPSE to perform isophotal model fits 
to the surface brightness distribution of the galaxy (after masking 
out all pixels covered by dust features and bright sources). Subtrac- 
tion of that model revealed a rich system of extensive shells and rip- 
ples (cf. also Schweizer 1980), as well as the dust features, small 
galaxies and compact sources. To also remove the shells and rip- 
ples, we first subtracted all compact sources using the DAOPHOT-II 
PSF fitting task. Then we applied a median filter (kernel = 9x9 
pixels) to the image. Finally, the sum of the isophotal model and 
the median-filtered image constituted the final model image. 

The sky background levels and their uncertainties were de- 
termined by fitting power laws to the outer parts of the intensity 
profiles of the isophotal model as described in Goudfrooij et al. 
(1994a). The selection and photometry of star cluster candidates 



was carried out using DAOPHOT-IL After aligning the B, V and 
Ic images on the same coordinate system, the DAOFIND task was 
run on the V image (again, divided by the square root of the model 
image) which had the highest signal-to-noise. The detection thresh- 
old was set at 4 sigma above the (residual) background. Due to the 
very complex dusty structure of the central regions of NGC 1316, 
especially at ground-based resolution, we decided to exclude the 
compact sources found within the central ~ 50" from the nucleus 
from further analysis. The great majority of those compact sources 
was measured on the WFPC2 images anyway. 

Aperture photometry of the cluster candidates was performed 
using an aperture radius of 4 pixels. Aperture corrections were de- 
rived for each filter from a curve-of-growth analysis of 10 well- 
exposed, isolated stars in the field. 

A _B vs. B — I CMD for the 350 detected objects is shown in 
Fig. dI The direction of the reddening vector is also shown. From 
this list of detections, a strict cluster candidate list was made by 
applying a colour selection criterion of 1.0 < B—I < 2.5, i.e., the 
same as for the WF frames of the HST data. Fig. ra is further dis- 
Photometry and astrometry for the 50 bright- 



cussed in Section 4. 1 



A2 in 



est cluster candidates in the NTT frames is shown in Table 

the Appendix. A similar table for the full list of cluster candidates 

in the NTT frames can be obtained in electronic form from CDS. 

Completeness tests were done the same way as for the 
HSTfWVFCl images (cf. Section 3.1). The resulting completeness 
functions for the B and V images are shown in Fig. n. The 50 per 
cent completeness levels are at _B ~ 24.10 and V ~ 23.85. These 
represent our faint magnitude cut-offs for the cluster candidates. 
Even using these selection criteria we do expect some contami- 
nation by foreground stars and background galaxies in the cluster 
candidate sample, due to the lower spatial resolution with respect to 



the HST data. This issue will be discussed below in Section 14 . 1 



3.3 Ground-based near-infrared imaging 

3.3.1 Observations and data reduction 

Near-infrared images of NGC 1316 were obtained during the nights 
of 1996 November 22-25 using the IRAC2 camera mounted on 
the ESO/MPI 2.2-m telescope. The detector was a Rockwell NIC- 
M0S3 array (HgCdTe, 256 x 256 pixels). We used lens C, having 
an image scale of 0'.'509 per pixel. The read-out noise was 22.5 e~ 
per pixel. We mapped a region of ~ 19 arcmin^ around the cen- 
tre of NGC 1316 with the J, H and K' filters. The field cover- 
ing the central region of NGC 1316 was exposed for a longer time 
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Figure 5. B vs. B~ I colour-magnitude diagram for NGC 1316 GC can- 
didates from the NTT data. Only objects with B — I errors less than 1 .0 
mag are plotted. The dotted lines delineate the colour and magnitude crite- 
ria for cluster candidates, while the hatched region represents the apparent 
magnitude and colour range for GCs in the halo of our Galaxy, placed at 
the distance of NGC 1316 (see text for details). The solid hne shows the 
reddening vector for Ag = 1 mag of extinction. 



than the other fields in order to get the highest possible signal-to- 
noise ratio in the region overlapping the WFPC2 field. The sky was 
photometric during the whole observing run, and the seeing varied 
between (X'S and I'.'O. Separate sky frames located at ^ 20' from 
the galaxy centre were also taken, interleaved with the individual 
galaxy frames to enable proper sky subtraction. The exposure times 
for each individual frame were chosen to be as long as possible, but 
care was taken to have the intensity at the galaxy centre stay within 
the 1 per cent linearity regime. 

The reduction procedure was as follows. First, all dark frames 
with identical integration time were combined and subtracted from 
the raw images (i.e., dome flats, standard star frames, sky frames, 
and galaxy frames). For each filter, the dome flats taken with iden- 
tical lamp voltages were combined, and the combined zero- voltage 
dome flat was subtracted from all other dome flats. We used a com- 
bined dome flat with an intensity level similar to that of the object 
frames for final flat fielding, after having normalized it to a mean 
value of unity. The dome flats were also used to create a bad pixel 
mask which allowed the replacement of bad pixels by the median 
of the values of the surrounding pixels. A final correction to the 
low-frequency flat field correction ('illumination correction') was 
made by mapping a standard star on a grid of 16 positions across 
the surface of the array, and fitting a third-order polynomial to the 
total intensities of the star at the different positions. 

The photometric calibration was derived from observations of 
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Figure 6. B and V band completeness functions for the NTT photometry 
of NGC 1316 cluster candidates. The 50 per cent completeness levels are at 
S~ 24.10 and V ~ 23.85. 

several standard stars from Carter & Meadows ( |199^ ). Each stan- 
dard star was observed on 2 regions of the array that are free of 
bad pixels. Aperture photometry was then performed on an image 
that is the difference of the two flat-fielded standard star images, 
ensuring a proper sky subtraction. The resulting colour equations 
were: 

J-j = (21.930 ±0.008) 

+ (0.148 ± 0.014) ij-k') - 0.15 x 
H-h = (21.706 ±0.006) 

± (0.028 ± 0.010) {j-k') - 0.10 a; 
K-k' = (21.170 ±0.006) 

± (0.022 ± 0.017) (j-fc') 



0.17a; 



where upper case JHK denote calibrated magnitudes, lower case 
jhk' denote instrumental magnitudes normalized to an exposure 
time of 1 s, and x is the airmass. These zeropoints were dou- 
ble checked by com paring the JHK aperture photometry of 
NGC 1316 by Glass ( 1984 ) with our data. They turned out to be 
consistent with one another to within 0.03 mag in any passband. 



3.3.2 Cluster candidate selection and photometry 

In order to put the optical and near-IR images on a common coordi- 
nate system, we measured the centroids of several point sources in 
common between the near-IR images and the NTT V-band image, 
using the IRAF task XYXYMATCH. The alignment error was always 
(and, mostly, far) below 0.1 pixel. All near-IR images were sub- 
sequently individually transformed to the NTT coordinate system. 
Finally, all transformed near-IR images were mosaiced together per 
filter. Fig. W shows the final _ff-band mosaiced image, with the loca- 
tion of the WFPC2 field superposed onto it. The galaxy background 
subtraction, point-source photometry and the derivation of aperture 
corrections were performed in the same way as for the optical NTT 
images. The near-IR photometry of the cluster candidates are listed 
in Tables Al and Wi in Appendix A. 



PROPERTIES OF THE STAR CLUSTER SYSTEM OF 
NGC 1316 
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Figure 7. Grey-scale plot of H-hand mosaic of the central 3.'75 X 3.'75 of NGC 1316, observed with IRAC2B at the ESO/MPI 2.2-m telescope (see Sec- 
tion B.3h. The location of the WFPC2 frames is outlined in black. The plate scale is indicated by the bar in the lower left of the image. North is up and East is 
to the left. The faint, extended 'object' located ~ 1 arcmin west of NGC 1316 is a detector feature (caused by amplifier crosstalk). 



4.1 Optical properties of stellar populations 

4.1.1 Interpretation of optical photometry 

The B vs. B - J CMD for cluster candidates in the WFPC2 and 
NTT frames was shown in Figs. W and H, respectively. The area 
which the GC system of the halo of the Milky Way would occupy 
if placed at the distance of NGC 1316 is also indicated, using the 
database of Harris (1996). These diagrams suggest some interesting 
results. 

First, the B—I colours of the clusters that are brighter than the 
brightest GC in our Galaxy are quite uniform and on the red end 
of the range covered by the Galactic GCs (and beyond). Spectra of 
three of these bright clusters were analyzed in Paper I, where we 
found all three to be ^ coeval with an age of 3.0 ± 0.5 Gyr and a 
solar metallicity (to within 0.15 dex). The B — I colours of these 
bright clusters are consistent with their spectroscopically derived 
age and metallicity range (Paper I; see also below). 

Second, the clusters fainter than B ~ 23 tend to be some- 
what bluer on average than the brighter clusters (excluding some 
inner clusters in the PC frame which are likely reddened by dust, 
cf. Fig. Bb. This trend is more obvious from the WFPC2 diagram 
(Fig. H) than from the NTT diagram (Fig. m, which is most likely 
due to the higher photometric accuracy of the WFPC2 data rela- 
tive to that of the NTT data. The range of colours encompassed 
by the clusters with B ^ 23 is roughly consistent with that en- 
compassed by the Galactic GCs, considering the photometric errors 
(e.g., aB-i ~ 0. 13 at B ~ 25). The main cause of the bluer mean 
colours at fainter brightness levels would seem to be a lower metal- 



licity rather than a younger age, as younger populations are brighter 
than older ones. These fainter clusters may be expected to be a mix- 
ture of (i) old, metal-poor clusters associated with the pre-merger 
galaxies (given their overlapping colours and luminosities) and (ii) 
the lower-luminosity end of the system of second-generation clus- 
ters formed during the merger. 

These results are illustrated further in Fig. plwhich shows B— / 
colour distributions of cluster candidates on the WFPC2 frames (in 
the form of histograms), and compares them with that of the GC 
system of the Milky Way. The colour distribution of the full cluster 
sample (depicted as open histogram on panels (a) and (b) of Fig. M) 
exhibits two peaks: one at B — J ~ 1.5 [coinciding with the peak 
of the Milky Way GC system, cf. panel (c)] and one at B — / ~ 1.8. 
Panel (a) shows clearly that the bright end of the cluster luminosity 
function (hereafter LF; shown as filled histogram; many of these 
objects were proven to be genuine clusters in Paper I) primarily 
makes up the red peak of the overall colour distribution, while the 
panel (b) shows that at the faint end of the cluster LF, the blue 
peak of the colour distribution becomes about as highly populated 
as the red peak. The reddest 'wing' to the colour distribution in 
Fig. H is mostly populated by cluster candidates in the innermost 
regions (i.e., the PC frame, cf. Fig. S). This suggests that the most 
likely main cause of the red wing is dust extinction rather than high 
metallicity. To decide unambiguously between these possibilities, 
spectroscopy and/or near-IR photometry with 8-lOm-class tele- 
scopes will be needed. Finally, the colour distribution of the clus- 
ter candidates on the NTT frames (outside a radius of 50", like in 
Fig. pi) is shown in panel (d) of Fig. M. It is similar to that of the 
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Figure 8. Histograms of the B — I colours of NGC 1316 cluster candidates 
compared with those of GCs in the Milky Way. In panels (a) and (b), the 
open histogram depicts the colour distribution of all cluster candidates in 
the WFPC2 field of NGC 1316. The shaded histogram in panel (a) depicts 
the colour distribution of the NGC 1316 cluster candidates with B < 24.0 
in the WFPC2 photometry, while the shaded histogram in panel (b) depicts 
the colour distribution of the NGC 1316 cluster candidates with B > 24.5 
in the WFPC2 photometry. Panel (c) depicts the histogram of B—I colours 
of GCs in the Milky Way. Panel (d) depicts the histogram oi B — I colours 
of NGC 1316 cluster candidates from the NTT photometry. 



clusters on the WFPC2 frames, except that the blue peak is less 
populated in the former. This is entirely consistent with our inter- 
pretation involving two cluster populations mentioned above, given 
the fact that the NTT photometry doesn't reach deep enough to de- 
tect the bulk of the fainter, bluer clusters that are supposedly the 
'old' clusters from the pre-merger galaxies. 

Summarizing the optical photometry results presented above, 
the blue peak of the overall B — I histogram, which coincides with 



the modus of the B — I colour distribution of the Galactic GC sys- 
tem, is primarily populated by the faint end of the cluster LF, while 
the red peak in the B—I histogram is solely populated by the bright 
end of the cluster LF. Since we also know from spectroscopy that 
the brightest clusters in NGC 1316 are 3 ± 0.5 Gyr old and have 
a solar metallicity to within 0.15 dex (Paper I), an obvious inter- 
pretation of these results is that we are dealing with the presence 
of two distinct cluster populations: (i) a second-generation popula- 
tion clusters of roughly solar metallicity that were formed ~ 3 Gyr 
ago during a major spiral-spiral merger (as determined from spec- 
tra in Paper I), and (ii) a population of old, primarily metal-poor 
clusters that were associated with the progenitor galaxies. From the 
colour histograms of the WFPC2 data shown in panels (a) and (b) 
of Fig. B, we estimate that of order 50 per cent of the total number 
of clusters is of the second generation. 

Is this interpretation consistent with stellar population model 
predictions ? We considered the following to address this ques- 
tion. The metallicity range encompassed by the brightest clus- 
ters at an age of 3 ± 0.5 Gyr (derived from spectra in Paper I: 
-0.15 <, [Z/Zq] <, 0.15) corresponds to 1.75 < B-I < 1.95 
according to the BC96 models. Considering the typical photomet- 
ric errors of the WFPC2 photometry as well (e.g., A (B — I) ^ 0.1 
for B < 24), we assigned a B-I interval of 1.65 < B-l < 2.05 
to the putative second-generation population of clusters. This is in- 
dicated as the hatched region on the left side of panel (a) of Fig. H, 
which depicts the predicted time evolution of the B — I colour for 
single-burst populations older than 1 Gyr [using the BC96 models 
as well as the Maraston (2001; see Sect. 4.2) models for compari- 



son]. The hatched region on the right side of panel (a) indicates the 
B — I interval occupied by the Galactic GC system which corre- 
sponds to a metallicity range of —1.5 ^ [Z/Zq] jS —0.5 (Harris 
1996; Minniti 1995) which is consistent with the model predictions 
at an age of 14 Gyr. These two B—I intervals are also indicated on 
panel (c) of Fig. |, which is the B vs. B-I CMD from the WFPC2 
photometry (i.e., the same as Fig.H). It can be seen that these two 
B—I intervals provide a very good fit to the positions of the cluster 
candidates on the CMD. The upper brightness limit to the hatched 
region for the second-generation population of clusters on panel (c) 
was chosen to encompass the brightest cluster candidates, since the 
spectra in Paper I showed them to be genuine clusters associated 
with NGC 1316. Note that this brightness limit is several magni- 
tudes brighter than that of the Galactic cluster system, rendering 
the brightest second-generation clusters in NGC 1316 to be much 
more massive than cu Cen, the brightest cluster in our Galaxy (this 
was discussed in detail in Paper I). 



4.1.2 Evolution of B — I colour distribution 

Will the colour distribution of the NGC 1316 cluster system evolve 
into a bimodal one, as typically found in giant ellipticals of similar 
luminosity? This is a key question in the context of finding evidence 
for or against the 'merger scenario' (Ashman & Zepf 1992) on the 
nature of bimodal colour distributions of GC systems of giant el- 
liptical galaxies. NGC 1316 constitutes an important probe in this 
context as well, since (/) the age and the metallicity of the bright- 
est second-generation GCs are known to relatively high accuracy 
(Paper I), and (//) the first signs of colour bimodality have already 
become apparent, with properties that are entirely consistent with 
the predictions of the merger scenario. 

So what will the colour distribution look like at a time when 
the signs of merger activity have faded away, e.g., after ^ 10 Gyr? 
This issue is addressed by panel (b) of Fig. M which depicts the time 
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Figure 9. Panel (a): Time evolution of B — I colour index according to single-burst stellar population models for ages old er th an 1 Gyr. The black lines 
represent the Brazual & Chariot (1996) models, while the grey lines represent the Maraston (2001) models (described in Sect. ^3). Model curves are plotted 
for a Salpeter (1955) IMF and the following metallicities: 0.02 solar (short-dashed lines), 0.2 solar (long-dashed lines), 1.0 solar (solid lines), and 2.5 solar 
(dot-long-dashed lines). Ages of 3 Gyr and 14 Gyr are indicated as vertical short-long-dashed lines. The region hatched by dotted lines on the left indicates 
the B — I colour interval populated by 3 Gyr old clusters in NGC 1316 with a solar metallicity to within 0.15 dex, while the hatched region on the right does 
so for a old, metal-poor cluster population such as that of the Milky Way halo [the same B — I intervals are indicated in Panel (c)]. 

Pane] (b): Evolution of B-band magnitudes (labels are absolute magnitudes on the left-hand side, apparent magnitudes at the distance of NGC 1316 on the 
right-hand side) for a globular cluster having the mass of oj Cen (5 X 10^ Mq, Meylan et al. 1995). Curves as in Panel (a). Ages of 3 Gyr. 10 Gyr. and 14 Gyr 
are indicated. The horizontal dotted lines indicate B magnitudes for specific cases, as discussed in Section [4. l| .2. 

Pane] (c): Same B vs. B—I colour-magnitude diagram as in Fig. hi The hatched region on the left (i.e., 1 .35 < B—I < 1 .75) delineates the apparent magnitude 
and colour range for GCs in the halo of our Galaxy (placed at the distance of NGC 1316), while the hatched region on the right (i.e., 1.65 < i? — J < 2.05) 
does so for a 3 Gyr old population of star clusters in NGC 1316 with a metallicity of [Z/Zq} = 0.00 ±0.15 dex. See discussion in Section Wjll. 



evolution of the absolute magnitude Mb (as well as B at the dis- 
tance of NGC 1316) for a globular cluster with the mass of uj Cen 
{M = 5 X 10'' M©, cf Meylan et al. 1995). One can see that at 
solar metallicity and an age of 3 Gyr, such a cluster is predicted to 
have B = 21.80, whereas it will fade to B = 23.05 at an age of 10 
Gyr, a fading of 1.25 mag. In that time span, its B — I colour will 
have reddened from ~ 1.85 to ~2.15 [cf. panel (a)]. The peaks in 
the B — I colour distribution at an age of 10 Gyr are thus expected 
to be at B — J = 1.55 and B — / = 2.15. This is remarkably similar 
to the peaks absented in bimodal distributions of cluster colours 
in well-studied giant ellipticals with luminosities similar to that of 
NGC 1316 (e.g., NGC 1399 and NGC 1404 (Forbes et al. 1998), 
NGC 3923 (Zepf, Ashman & Geisler 1995b), NGC 4472 (Geisler, 
Lee & Kim 1996), NGC 4486 (Kundu et al. 1999). Hence, the an- 
swer to the question posed as the first sentence of this subsection is 
affirmative. The cluster system of this merger remnant will evolve 
into one with a typical bimodal colour distribution (barring further 
merger events), providing strong support in favour of the 'merger 
scenario' . 



4.1.3 Radial distribution of the two subpopulations 

Another key prediction of the 'merger scenario' is that the newly 
formed, metal-rich clusters will be more centrally concentrated 
than the pre-existing metal-poor clusters. However, the multi -phase 
collapse (non-merger) scenario proposed by Forbes et al. (1997) 



predicts a very similar distribution of the two subpopulations, so 
that the occurrence of this property in 'normal' giant ellipticals 
featuring bimodal colour distributions does not provide strong ev- 
idence in support of either scenario. As discussed in the previous 
subsection however, GC systems of intermediate-age merger rem- 
nants such as NGC 1316 are important probes in this context. We 
compare the radial distribution of the two subpopulations by using 
the WFPC2 data, assigning clusters with 1.00 < B- J < 1.55 and 
A(B — /) < 0.1 to represent the 'metal-poor' ones and clusters 
with 1.75 < B-I < 2.5 and A{B - 1) < 0.1 to represent the 
second-generation 'metal-rich' ones. These criteria ensure that the 
B — I interval that overlaps with both subpopulations is excluded 
from this exercise. The radial distribution of both subpopulations is 
depicted in Fig. hol A glance at the lower panel (a diagram of B—I 
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Figure 10. Radial distribution of the two cluster subpopulations from the 
WFPC2 data. The lower panel shows B — / as a function of galactocen- 
tric radius, while the upper panel shows surface density as a function of 
galactocentric radius. The surface densities were measured by averaging 
over three logarithmically spaced annuli. Filled squares represent the pro- 
posed 'metal-rich', second-generation clusters and open squares represent 
the proposed 'metal-poor', pre-merger generation of clusters. Error bars in 
the lower panel are only plotted if larger than the symbol size. The surface 
densities have been c oiTe cted for foreground stars and background galaxies 
as discussed in Sect. |t.4ll. 



vs. radius), already shows clearly that the 'metal-rich' clusters are 
more centrally concentrated than the 'metal-poor' ones. This con- 
clusion is endorsed by the upper panel which shows a plot of clus- 
ter surface density vs. radius, derived after counting the number of 
clusters within three radial annuli, logarithmically spaced in galac- 
tocentric radius. Note that a similar result was found for NGC 3610, 
a candidate 4 Gyr old merger remnant featuring a similarly emerg- 
ing colour bimodality (Whitmore et al. 1997). 



4.2 Near-IR properties of stellar populations 

In interpreting the JHK photometry of star cluster candidates in 
NGC 1316, we restrict ourselves to the objects that have optical 
photometry as well. The near-IR photometry was performed in part 
because it serves as a very sensitive check for the presence of stel- 
lar populations younger than ~ 2 Gyr, which is relevant to the is- 
sues at hand in this paper. Near-IR photometry is remarkably sen- 
sitive to the presence of thermally pulsing asymptotic giant branch 
(TP-AGB) stars. TP-AGB stars are known to undergo a rather ac- 
tive phase during which a large amount of fuel in extended hy- 
drogen and helium shells around degenerate carbon-oxygen cores 
is burned, which results in a sharp increase of the near-IR stellar 
luminosity, particularly in the A^-ban d. This so-called 'AGB phase 
transition' ( Renzini & Buzzoni 198f ) starts at ages ~ 200-300 Myr 
and lasts > 



1 Gyr (cf. Renzini & Voli 1981; Iben & Renzini 1983; 
Maraston 1998; Langon 1999). 

In contrast to near-IR photometry, the evolution of optical 
colour indices such as B — 1^ or B — J for a simple stellar popula- 
tion is monotonic with age, since they mainly trace the evolution of 
the turn-off part of the population. Hence, colour-colour diagrams 
involving such an optical colour index along with an index that 
is sensitive to the AGB phase transition (e.g., V — K, I — K, or 



J — K) can be expected to provide powerful means of identifying 
intermediate-age stellar populations. 

However, the actual inclusion of the TP-AGB contribution in 
population synthesis models is complicated by poor knowledge of 
the physical mechanisms driving AGB evolution. Specifically, the 
amount of fuel burned during the thermally pulsing phase is largely 
unknown due to uncertainties in the mass loss rate, mixing, and 
the efficiency of hydrogen burning at the bottom of the convective 
layer (Renzini & Buzzoni 1986). Fortunately, intermediate-age star 
clusters in the Magellanic Clouds cover the relevant age range (see, 
e.g., the review by Olszewski 1995) and thus offer an oppor tunity 
to calibrate these effects. This was done by Maraston (1998) who 
presented SSP models in which the contribution of TP-AGB stars to 
the total SSP energy as a function of age is calibrated on observed 
photometry of star clusters in the Magellanic Clouds. Maraston et 
al. (2001) computed similar models for va rious metallicities, adopt- 
ing the prescriptions by Renzini & Voli (1981) to describe the in- 
fluence of metallicity on the TP-AGB energy. These models are 
described in the next subsection (see Maraston et al. 2001 for more 
details). 

4.2.1 SSP modelling of populations including TP-AGB stars 

We used the recent sets of SSP models with metallicities of 0.5 
Zq, 1.0 Zq, a nd 2.0 Zq that are described in detail in Maras- 



ton et al. ( 2001 ) and Maraston (2001, in preparation). In summary. 



these models were c omputed with the evolutionary synthesis code 
of Maraston (1998), which is based upon the fuel consumption 



theorem (Renzini & Buzzoni 1986). The AGB contribution to the 
bolometric luminosity as a function of age was calibrated using 
the observed AGB contributions in Magellanic Cloud clusters from 
Frogel, Mould & Blanco (199G), grouped into age bins adopting 



the classification of Searle, Wilkinson & Bagnuolo (1980, hereafter 
SWB). Model colours were computed by calibrating the partition 
of the total TP-AGB fuel among stars of spectral type C (Carbon 
stars) and M. The fractional number of C stars relative to the total 
number of AGB stars, Nc, is a function of SSP age and metallicity. 
The age dependence o f Nc was evaluate d from the data for Mag- 
ellanic Cloud clusters ( Frogel etal. 199C| ), while the metallicity de- 
pend ence o f Nc was evaluated using the prescription of Renzini & 
Voli (1981 ). Observed co lours for C- an d M-t ype AGB stars were 



tak en fro m Cohen et al. ( |1981[ ), Cohen ( |1982| ), and Westerlund et 
al. (|l99l|). 



The success of these new models in describing the photo- 
metric properties of intermediate-age star clusters was shown in 
Maraston et al. (2001, especially their Fig. 5). The star clusters in 
NGC 7252, a merger remnant galaxy of age ~ 300-500 Myr (e.g., 
Hibbard & Mihos 1995; Schweizer & Seitzer 1998), were found 
to exhibit very red V ~ K colours for their B — V, which is ex- 
actly what is observed for Magellanic Cloud clusters of SWB types 
IV-VI. Maraston et al. (2001) showed clearly that this behaviour 
cannot be understood without the addition of the TP-AGB contri- 
bution, while it was readily fit by the new models for 0.5 Zq. 

Given the complex and still poorly determined shape of the 
spectrum of TP-AGB stars, intermediate-age SSP colours must be 
calibrated against observations in every band. For the Cousins I- 
band, observations of C- and M-type TP-AGB stars are not avail- 
able, nor are integrated / magnitudes of globular clusters in the 
Magellanic Clouds. For the purpose of the present paper, the V — I 
model colours have therefore been calibrated on observations of 
intermediate-age star clusters in NGC 7252 (from Whitmore et al. 
1997; see Maraston et al. 2001, Table 2). The dereddened V-I 
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Figure 11. Optical-IR colour-magnitude diagrams of cluster candidates [B vs. B—K in (a); K vs. B—K in (b)]. Symbol types as in Fig. BJ The symbol sizes 
have been chosen to scale linearly with projected distance to the centre of NGC 1316. Error bars are only plotted if they are larger than the symbol size. The 
grey-scale region at the bottom left delineates the apparent magnitude and colour range for GCs in our Galaxy, placed at the distance of NGC 1316. The solid 
arrow indicates the reddening vector for Ag = 1 mag of extinction, while the dashed arrow indicates the vector for a metallicity increase of 0.33 dex at an age 
of 3 Gyr (see text for details). 



colours of these clusters span a small range: V — I = 0.64 ± 0.04. 
Their ages have been found to be ~ 300 Myr (Maraston et al. 2001, 
consistent with the results by Schweizer & Seitzer 1998 based on 
optical spectroscopy). Therefore the model V — 7 at age 300 Myr 
has been calibrated to reproduce V — I = 0.64. Note that this pro- 
cedure retains the amount of fuel that was calibrated in Maraston 
(1998). The calibrating variable has been chosen to be the V — I 
colour of the input TP-AGB stars. Specifically, V — I = 1.1 for TP- 
AGB M stars and 1.4 ^ 1/-/ ^ 2.0 for C-type TP-AGB stars is 
found to perfectly reproduce an integrated colour V — I = 0.64. 

These same input colours for TP-AGB stars have been used 
to compute the models older than 300 Myr. This is a reasonable 
assumption because the colours of C stars are not found to change 
significantly with the age of the parent star cluster (Frogel et al. 
1990 and references therein), the evolving parameter being the fuel 
consumption. 



4.2.2 Results 

Fig. [n| shows the B vs. B - K and K vs. B - K CMDs for the 
detected cluster candidates in NGC 1316. Their location in these 
diagrams is compared to the area covered by the system of halo 
GCs in our Galaxy (shown in grey scale). The JHK photometry 
of the Galactic GCs was taken from Aaronson et al. (1978, and 
private communication). Only data with photometric errors in B—K 
smaller than 1 mag are plotted. The detected star clusters are more 
luminous (as well as redder) than the halo GCs in our Galaxy, and 
are thus likely belonging to the second-generation cluster system 
formed during the merger. Comparing the two CMDs in Fig. hlL 
one can see that the colour of the reddest inner clusters (the open 
circles; these clusters are located within ~ 20" from the nucleus) is 
most likely due to extinction within the dusty filaments rather than 
being due to higher metallicity. The present near-IR imaging does 
not reach deep enough to detect the clusters belonging to the 'blue 
peak' in the B—I histogram (cf. Figs. |8|and|9|). 8-m class telescopes 
are needed for this purpose, which is important to obtain highly 
increased power (with respect to optical colours) in disentangling 
different cluster populations in galaxies. 



Fig. Il3 shows the B~I vs. J—K colour-colour diagram of the 
cluster candidates in NGC 1316. This diagram is particularly sensi- 
tive to the presence of TP-AGB stars, as can be evaluated from the 
models of Maraston (2001, described above) which are overplotted 
in the left panel for metallicities of 0.5 Zq, 1.0 Zq, and 2.0 Zq. 
Symbols along the model curves indicate a number of reference 
ages. For comparison, the right panel shows the same data along 
with curves describing the BC96 models (which do not include the 
calibrated effect of the AGB phase transition). Fig. O shows clearly 
that the colours of the best-observed clusters in NGC 1316 are in- 
deed not well described by an intermediate-age population domi- 
nated by TP-AGB stars such as the clusters of NGC 7252. Instead, 
the NGC 1316 clusters are best described by the solar metallicity 
model at an age of 2.5-3 Gyr (allowing for reddening of the red- 
dest inner clusters by dust as seen before), again consistent with 
the spectroscopic results of Paper I and the optical photometry dis- 



cussed in Section %. 1 



4.3 Star cluster luminosity function 

The typical luminosity function (LF) of 'old' globular clusters 
around galaxies such as our Milky Way or 'normar giant ellipti- 
cal galaxies such as NGC 1399 or M87 is well represented by a 
log-normal function (i.e., a Gaussian in magnitude units; see e.g.. 
Ashman & Zepf 1998 and references therein) of which the peak 
absolute magnitude concides with that of the GC system of our 
Galaxy (e.g., Whitmore 1997). In contrast, the shape of luminosity 
functions of clusters in 'young' [age ^ 500 Myr] merger remnant 
or starburst galaxies is a power law with a slope of ~ — 2 (see, e.g., 
Meurer et al. 1995; Schweizer et al. 1996; Miller et al. 1997; Whit- 
more et al. 1999b; Zepf et al. 1999). This decidedly different be- 
haviour has been used to argue that the nature of young clusters in 
mergers or starburst galaxies is fundamen tally different from clus- 
ters in old systems (van den Bergh 1995). On the other hand how- 



ever, several theoretical investigations have shown that dynamical 
evolution can significantly alter the mass function of star clusters 
over time (e.g.. Fall & Rees 1977; Aguilar, Hut & Ostriker 1988; 
Meylan &. Heggie 1997; Murali & Weinberg 1997; Gnedin & Os- 
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Figure 12. B—I vs. J~K colour-colour diagram for star cluster candidates in NGC 1316. Symbol types and sizes of the cluster data points as in Fig.Mll. The 
grey-scale region at the bottom left of each panel delineates the observed colour ranges for GCs in our Galaxy. The solid arrow indicates the reddening vector 
for Ag = 1 mag of dust extinction. The left panel (a) shows the data along with curves and symbols depicting the Maraston (2001) models for a Salpeter 
(1955) IMF, for the metallicities 0.5 solar (dashed line and open triangles), 1.0 solar (solid line and filled squares) and 2.0 solar (dot-long-dashed line and filled 
triangles). Symbols along the model curves are plotted for ages of (from right to left) 14, 10, 4, 3, 2, 1, 0.8, and 0.6 Gyr The Maraston model for 0.02 solar 
metallicity and 15 Gyr age is shown by an open square for comparison. The right panel (b) shows the same data along with curves depicting the BC96 models 
for a Salpeter IMF. Model metallicities, line types, and symbols are the same as for panel (a). The BC96 model for 0.02 solar metallicity and 15 Gyr age is 
shown by an open square for comparison. For the BC96 models, some of the symbols corresponding to the youngest ages fall off the plot. 



triker 1997). In fact, dynamical evolution has been crucial in shap- 
ing the LFs of old cluster systems as they are today (cf. Gnedin & 
Ostriker 1997). In order to study the evolution of cluster LFs by 
observation as well as theory, it is important to find systems with 
ages that are intermediate between those of the young mergers and 
the ^ 10 Gyr old systems. NGC 1316 is an interesting probe in this 
respect, as the bright GCs in NGC 1316 are ~ 3 Gyr old (Paper I). 
For instance, if the GC system of NGC 1316 indeed is a combina- 
tion of "old" GCs that belonged to the merger progenitor galaxies 
and second-generation GCs that were formed in the merger, one 
would expect to find different LFs for the two subpopulations. 

Given the magnitude depth of the HSTfWFPCl data, its su- 
perb power in discarding background galaxies as star cluster can- 
didates, and the availability of a background field, we use the 
HST data for examining the star cluster LF for NGC 1316. While 
this it em was already discussed to some extent by Grillmair et al. 
(1999), we are interested in the exact functional representation of 
the shape of the LF for reasons outlined above. The LF for the 
cluster system of NGC 1316 is plotted in Fig. hsl Four panels 
are shown: One for all 'red' GCs, one for all 'blue' GCs (where 
'red' and 'blue' are defined as for Fig. hd in Sect. kj|.3), one for 
the full list of cluster candidates found on all WFPC2 chips, and 
one for cluster candidates located at radii > 50" from the centre 
(i.e., outside the dusty central filaments; hereafter called the 'un- 
obscured sample'). The solid histograms depict the observed LFs, 
while the dashed histograms depict the completeness-corrected and 
background- subtracted LFs. The completeness-corrected LF of the 
full sample of GCs in NGC 1316 can be seen to continue to increase 
out to beyond the 50 per cent completeness limit of the WFPC2 
data (cf. also Grillmair et al. 1999), reminiscent of the LFs found 
in young merger remnant galaxies as mentioned above. However, 
the two upper panels in Fig. lo clearly show that this behaviour is 
mostly due to the 'red' GCs, as the LF of the 'blue' GCs shows 
strong evidence for a 'turnover', similar to log-normal LFs such as 
those found for GCs in our Galaxy and in 'normal' ellipticals. For 
reference, the absolute magnitude of the 'turnover' of the Galac- 



tic GC system is Mb ~ —6.6 (Harris 1996), similar to the peak 
in the LF of the 'blue' GCs in NGC 1316. This is consistent with 
the notion that the 'blue' GCs in NGC 1316 stem from the merger 
progenitor spirals. Given the fact that the turnover is near the 50 
per cent completeness level however, this important result should 
be confirmed with deeper photometry. 

Given the well-known power-law nature of LFs of young and 
mixed-age cluster systems (see, e.g., Elson & Fall 1985 for the 
LMC; Whitmore et al. 1999b for NGC 4038/4039), we fitted the 
completeness-corrected, background-subtracted LFs for NGC 1316 
by a power law of the form 

0(L) dL (X L" dL. 

To limit the impact of uncertainties in the incompleteness factors, 
the fits were limited to bins containing objects with B < 25.6. 
Fig. |l3 shows the results of these fits (solid curves) performed by 
least squares with weights proportional to A'^ (the number of ob- 
jects per bin). The best fit for the 'red' GC candidate sample yields 
a = -1.23 ± 0.26. For the 'full' GC candidate the fit yields 
a = —0.97 ± 0.09, while that for the unobscured sample yields 
a = —1.30 ± 0.16. These LF power-law slopes are significantly 
flatter than those fou nd in 'young' mergers (cf. a — —2.6 ±0.2 
for NGC 40 38/4039 (|Whitmore et al. 1999b[ ), a = -2.1 ± 0.3 for 
NGC 3921 dSchweizer et al. 1996| ). a =^8±0.1 for NGC 7252 
( [Vliller et al. 1997]), and a ^ -1.8 for NGC 3256 (Z epf et al. 
1999)) and starburst galaxies {a ~ —2.0, cf. Meurer et al. 1995). 



4.3.1 Dynamical evolution of star cluster mass function 

The above results indicate that the LF (and mass function) of a star 
cluster system of a ~ 3 Gyr old merger remnant such as NGC 1316 
is a power law, with a slope that is significantly flatter than star 
cluster systems of young merger remnants. We suggest that this 
flattening of the power-law LF may represent an intermediate stage 
in the dynamical evolution of star cluster mass functions. We illus- 
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Figure 13. B-band luminosity functions of star cluster candidates using 
HST data. The top panel (a) shows the luminosity function for the red clus- 
ter candidates on the WFPC2 frames, panel (b) does so for the blue clus- 
ter candidates (see discussion in Section pt.3| on the definition of 'red' and 
'blue' clusters in this context), panel (c) does so for all cluster candidates 
together, and the bottom panel (d) does so for the star clusters found outside 
the dusty central filaments. Bottom axis unit is apparent magnitude B, and 
top axis unit is absolute magnitude Mg at the distance of NGC 1316. The 
histogram is filled for magnitudes below the 50 per cent completeness limit, 
and open beyond it. Dashed lines show the effect of dividing by the com- 
pleteness fractions and subtracting the (completeness-corrected) luminosity 
function of the background field (cf. Elson et al. 1998). Solid smooth curves 
are power-law fits to the completeness-corrected luminosity functions. 



trate this below using simple scaling relations between time scales 
of known long-term cluster disruption processes. 

The two dynamical processes that are most relevant in deter- 
mining cluster mass functions are (i) (internal) evaporation through 
escaping stars and (ii) tidal shocking. The characteristic time scales 
for these destructive processes are: 



0.138 
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Gi/2 m. ln(A) 



(1) 



for evaporation (Spitzer & Hart 1971; Gnedin & Ostriker 1997), 
and 
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for tidal shocking (Aguilar et al. 1988). In Eqs. (1) and (2), Mc is 
the total cluster mass, Vh is the cluster half-mass radius, G is the 
gravitational constant, m, is the average stellar mass in the cluster, 
ln(A) = In (0.4A'^) is the Coulomb logarithm (A'^ being the num- 
ber of stars in the cluster), Pr is the radial period of the cluster 
orbit, -Rpcri is the perigalacticon distance, Vc is the cluster velocity 
at perigalacticon, A^g is the galaxy mass, x(^peri) is a dimen- 
sionless factor that takes the mass distribution of the galaxy into 
account (cf. Aguilar & White 1985), (f^ ) is the mean square clus- 
ter radius, 77* is the net shock efficiency, ^ is the ratio of the frac- 
tional change of cluster binding energy to cluster mass induced by 
the shock (Spitzer & Chevalier 1972), and A(e, -Rperi) is a dimen- 
sionless correction factor proportional to the square of the differ- 
ence in force at perigalacticon vs. at apogalacticon (cf. Eq. (11) in 
Aguilar et al. 1988). For the purpose of this comparative exercise, 
we neglect the dependence of Tshock on the cluster orbit and galaxy 
potential, i.e., we assume that we are comparing cluster systems in 
(giant) ellipticals covering a relatively small range of galaxy masses 



(and cluster orbits). This leaves us with Tc- 



« M 



1/2 3/2 



and 



Tshock CC McT-^ . 

Imaging studies of young clusters in merger remnants and 
starburst galaxies using HST have revealed a very weak dependence 
of cluster radii on luminosity (e.g., r^ oc l,~0'^^ in NGC 3256, 
Zepf et al. 1999; cf. also Meurer et al. 1995). Using Zepf et al.'s 
results, the relative scalings of destruction time scales as a function 
of cluster mass become Tcvap oc Ale'® and Tshock oc AlJ,'*. These 
relations can be used to determine the time evolution of any charac- 
teristic cluster mass that is undergoing a particular stage of destruc- 
tion like e.g., the turnover mass of 'old' cluster systems in galaxies 
such as the Milky Way or nearby 'normal' elliptical galaxies. As- 
suming an age of 14 Gyr for the old cluster systems, the turnover 
mass at an age of 3 Gyr (i.e., the age of the second-generation clus- 
ters in NGC 1 3 16) will be 8 per cent of that in the Milky Way if 
evaporation is the dominant destruction mechanism, or 15 per cent 
of that in the Milky Way if tidal shocking is the dominant destruc- 
tion mechanism. A relevant question to ask in the context of the ob- 
served LF of clusters in NGC 1 3 16 is then: Would one be able to de- 
tect clusters with masses of 8 - 15 per cent of the turnover mass of 
the Milky Way GC system in the WFPC2 data? According to recent 
stellar population models and a Salpeter (1955) IMF, the M/ Lb 
ratio of a 3 Gyr old stellar population of solar metallicity is a fac- 
tor 1.61 lower than that of a 14 Gyr old population with 0.05 solar 
metallicity (Maraston 1999; the factor is 1.44 when using the BC96 
models)rl Given the turnover magnitude of Mb ~ —6.6 in the 



+ The dependence on IMF is negligible in this context: E.g., For a Scalo 
(1986) IMF, the factor is 1.57 rather than 1.61. 
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Milky Way, one would then expect to see the turnover for the 'red' 
GCs in NGC 1316 occur within the range -5.1 <, Mb ^ -4.4, 
where the precise absolute magnitude of the turnover mainly de- 
pends on the relative importance of evaporation vs. tidal shocking. 
Note that this is beyond the 50 per cent completeness limit of the 
existing WFPC2 photometry (cf. Fig. O). Further observations of 
NGC 1316 are necessary to test whether the turnover in the cluster 
LF indeed appears at the expected luminosity. Resolving this issue 
should be quite easy with the Advanced Camera for Surveys (ACS) 
which is currently scheduled to be commissioned in 2002 aboard 
HST. 



4.4 Radial distribution and specific frequency 

The specific frequency of GCs, 



(Harris & van den Bergh 1981 ), i.e., the number of star clusters per 
galaxy luminosity in units of My = —15, is an important param- 
eter of star cluster systems of galaxies. Sjv is known to increase 
systematically along the Hubble sequence, from (Sjv) =0.5 ± 0.2 
for Sc spirals to 2.6 ± 0.5 for ellipticals outside galaxy clus ters, 
albeit with considerable galaxy-to-galaxy scatter (Harrisl991). As 
mentioned in Sect. mI a sizable (and still steadily growing) number 
of HST observations of merger and starburst galaxies have led to 
a wealth of discoveries of young star clusters, consistent with the 
predictions of the 'merger scenario' for giant ellipticals (Schweizer 
1987; Ashman & Zepf 1998 and references therein). An interesting 
question is then: Can the difference in specific frequency between 
ellipticals and spirals be accounted for by secondary populations 
of clusters created during (gas-rich) mergers? Again, NGC 1316 
is a well-suited probe to test this possibility, in the sense that we 
know that its body has already settled down to a typical r^' * sur- 
face brightness law, we know the age of the major merger during 
which the secondary population of clusters were created, and it oc- 
curred long ago for dynamical evolution of the cluster system hav- 
ing already had a significant impact, making estimates for its future 
evolution less uncertain than for younger mergers. In the remainder 
of this section, we derive the current value of Sn for NGC 1316 
and place limits on its future evolution. 



4.4.1 Radial surface density distribution 

In order to determine the radial surface density distribution and the 
total number of clusters in NGC 1316, we use a hybrid method 
that combines the depth of the WFPC2 imaging with the much 
larger sky area covered by the NTT imaging. The resulting radial 
surface density distribution is shown in Fig. [14, and was derived 
as follows. We use the NTT imaging results for clusters outside 
the area covered by the dusty filaments (beyond a radius of 50"), 
while we use the WFPC2 results inside that radius. Only clusters 
brighter than the NTT incompleteness limit were used (B < 24.0). 
In both datasets, we calculated the number of clusters in four an- 
nuli, logarithmically spaced in galactocentric radius. Next, we di- 
vided the number of clusters by the appropriate area coverage to 
yield a surface density (taking into account the limited azimuthal 
sky coverage of the WFPC2 data). Finally, we have to subtract the 
surface density of foreground stars and compact background galax- 
ies. For the WFPC2 data, we used the detected point sources in 
the field published by Elson et al. (1998; see also Grillmair et al. 
1999) who observed a background field in the Fornax galaxy clus- 
ter with WFPC2, using the same filters and exposure times as in 



the NGC 1316 WFPC2 observations. For the NTT data, we con- 
sidered the surface density of point sources in the WFPC2 back- 
ground field to be a lower limit to the surface density of foreground 
stars and background galaxies in the NTT cluster candidate list, 
since ground-based seeing doesn't allow one to discriminate point 
sources from distant compact galaxies. We estimated the level of 
background galaxy contamination using the Canada-France Red- 
shift Survey (CFRS) data (Lilly et al. 1995), within the appropriate 
region of the B, B — I CMD. Background galaxies at the bright 
end of our B-magnitude range peak at B — / ~ 1 .0 and become 
bluer towards fainter magnitudes. The resulting estimate of the sur- 
face density of background galaxies down to B = 24.0 is 1.3 ± 0.6 
arcmin"^ (cf. also Elson et al. 1998). The background-corrected 
surface densities of star clusters from the NTT data (hereafter re- 
ferred to as 'case A') are shown as filled circles in Fig. h4t open 
circles depict the surface densities that are only corrected for the 
detected targets in the WFPC2 background field (hereafter referred 
to as 'case B'). Since some fraction of background galaxies in the 
NTT field has already been removed by means of our DAOPHOT-II 
sharpness and roundness selection criteria, the true surface densi- 
ties will lie somewhere in between the open and filled symbols. 

Fig. h4| depicts the surface density distribution of the clusters, 
as well as that of the integrated light of the galaxy for comparison. 
The cluster surface density is seen to flatten off (in log-log space) 
towards the centre, relative to the integrated light profile (note that 
this excludes the possibility of being due to the presence of dust 
features in the inner regions). This is a known feature of GC sys- 
tems of 'normal, old' giant ellipticals, in which this flattening is 
somewhat more pronounced (e.g., Grillmair et al. 1994; Forbes et 
al. 1998). This feature is most probably due to the accumulative 
effects of tidal shocking (which is most effective in the central re- 
gions, cf. Eq.pl) and dynamical friction during the first few dynami- 
cal time-scales after the assembly of the galaxy (see, e.g., Gnedin & 
Ostriker 1997). The somewhat more moderate flattening of the sur- 
face density of the cluster system of NGC 1316 towards the centre 
relative to that of 'old' giant ellipticals is consistent with the notion 
that NGC 1316 is only a few Gyr old. 

Excluding the data points within a galactocentric radius of 
40", we made least-square fits to the surface density data using 
a power law of the form 

p(r) = por°. 

The best fit for the background-corrected data ('case A') is 
parametrized by po = (13.0 ± 1.1) arcmin"^ and a = —1.61 ± 
0.08, and is depicted by a solid line in Fig. h4|. The corresponding 
best-fitted slope for the data that were only corrected for the counts 
in the WFPC2 background field ('case B') is q = -1.36 ± 0.10, 
whereas the slope for the integrated B-band light of the galaxy is 
Q = — 1.41 ± 0.02 in the same radial range. We conclude that the 
radial surface density of the clusters (down to B = 24.0) is con- 
sistent with that of the underlying galaxy light. This situation is 
somewhat different from the case of GC systems of 'normal' giant 
ellipticals, for which the cluster distribution is typically somewhat 
more extended than that of the underlying starlight (e.g., Fleming 
et al. 1995; Forbes et al. 1998; Kundu et al. 1999). This is once 
again consistent with the intermediate-age nature of the NGC 1316 
cluster system in view of the dynamical evolution-related issues 
mentioned above. 
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Figure 14. Radial surface density profiles for the star cluster system of 
NGC 1316 down io B = 24.0. Filled squares represent measurements from 
the WFPC2 data, corrected for foreground stars and background galaxies as 
explained in the text. Circles represent measurements from the NTT data; 
the open circles have been coirected for foreground stars and background 
galaxies in the same way as the WFPC2 data, while the filled circles have 
been corrected for the density of background galaxies as estimated from 
the CFRS catalog (see text for details). The dot-dashed line depicts this 
estimated density of background galaxies, and its uncertainty is indicated 
by the errorbar at the bottom right of the figure. A power-law fit to the 
background-corrected data beyond the inner 'core region' is shown by a 
solid Une. The surface brightness profile of the integrated B-band light of 
the galaxy (normalized to the surface density of star clusters at a radius of 
1 arcmin) is shown by the dashed line. 



4.4.2 Present-day specific frequency 

Using the radial surface density profile, we can estimate the total 
number of clusters A^'gc and the specific frequency Sn- In view of 
the observed radial extent of the surface density distribution, we de- 
cided to use 10 ± 1 arcmin as the outer radius for the cluster system 
(the uncertainty in the radial extension of the cluster system was 
chosen somewhat arbitrarily, but in accordance with numerous pre- 
vious studies of cluster systems). Given the power-law nature of the 
LF of the cluster system of NGC 1316, the derivation of Sjv is not 
as straightforward as for systems featuring log-normal LPs. In par- 
ticular, the derived Sn is a lower limit, as it depends on the applied 
faint-magnitude limit. Integrating the power-law fit to the surface 
density between log (radius) = —0.2 and -l-l.O ± 0.1 (with radius 
in units of arcmin), integrat ing o ver the completeness-corrected LF 
(described above in Section p3| ), and adding the number of clusters 
found in the WFPC2 data interior to log (radius) = —0.2 yields the 
total number of clusters. Table 111 lists the total number of clusters 
down to i? = 25.5 and B = 26.0. (beyond which the completeness 
corrections become too uncertain) both for cases A and B described 



Magnitude limit 



B < 25.5 
B < 26.0 



JVgc ■S'jv 

(Case A)* 



1910 ± 120 
2405 ±210 



1.2 ±0.2 
1.5 ±0.3 



A^GC -S'jv 
(Case B)* 



2805 ± 225 
3535 ± 400 



1.5 ±0.2 
1.9 ±0.3 



Note to Table m * Case A represents the case in which estimated contami- 
nation from both foreground stars and background galaxies (the latter esti- 
mate uses data from the CFRS survey) has been corrected for, while Case 
B is only connected for the estimated contamination from foreground sjars. 
Reality is likely to lie in between the two cases, as discussed in Sect. \^.^ 

above. The corresponding specific frequencies are also listed, as- 
suming VS = 8.53 (RC3), equivalent to Mv = -23.27. 

As discussed above, the 'real' specific frequency is most likely 
to lie in between the results for cases A and B. Hence, S'jv = 1.7 
± 0.4 is our best estimate for the present-day specific frequency 
of clusters with B < 26.0. While this value of Sn is on the low 
side for giant ellipticals, it is by no means extraordinary low: e.g., 
6 giant ellipticals in the compilation of Harris (1991) have values 
lower than 1.7. 



4.4.3 Evolution of the specific frequency 

In order to draw a parallel between the 5jv in NGC 1316 and that 
of 'normal, old' giant ellipticals (e.g., with a luminosity-weighted 
population age of ~ 10 Gyr), one has to take into account three 
main evolutionary effects: (i) Luminosity fading of the integrated 
V-band light of NGC 1316, (ii) luminosity fading of the 3-Gyr-old, 
second-generation clusters, and (/(/) disruption of low-mass clus- 
ters in the potential well of NGC 1316. As to effect (i) above, re- 
cent studies have shown that the (luminosity-weighted) age of the 
integrated light of NGC 1316 within ~ 20" radius is ~ 2 Gyr, as de- 
rived from optical line-strength indices (Kuntschner 2000). The age 
of the integrated light outside 20" radius is not as well constrained. 
However, the optical and near-IR colours in those regions are sim- 
ilar to those of 'normal' giant ellipticals (Schweizer 1980; Caon et 
al. 1994; Glass 1984), which may indicate that the stellar popula- 
tion there is similarly 'old'. Table H lists the evolution of Aly with 
time [using the Maraston (1998) SSP models for solar metallicity 
and a Salpeter IMF; the BC96 models give very similar results] for 
two simplified situations: (a) assuming that the current age of the 
integrated light of the whole galaxy is 2 Gyr, and (b) assuming that 
the age of the integrated light in the inner 20" is 2 Gyr, and 14 Gyr 
outside that. We consider reality to be most likely nearer the latter 
assumption. As to effect (ii) above, we considered the following. 
We divided the cluster candidate sample up into 'blue' and a 'red' 
subsamples. In this context, the 'blue' subsample was defined as 
the candidates with 1.0 < B — / < 1.65, whereas the 'red' subsam- 
ple was defined by 1.65 < B — 7 < 2.5. The B magnitudes of the 
'red' clusters were then faded to an age of 10 Gyr according to the 
Maraston (2001) SSP models for solar metallicity and a Salpeter 
IMF. A new LF was created for this evolved 'red' GC population, 
which is shown in Fig. flSl To estimate the expected turnover mag- 
nitude for this GC population, we considered that for a GC of a 
given mass, the B magnitude for solar metallicity is 1. 1 mag fainter 
than for 0.03 solar (the median metallicity of the GC system of our 
Galaxy). Hence, the turnover magnitude for a solar metallicity GC 
system should be at Mb ~ —6.6 -F 1.1 = —5.5 mag at 14 Gyr age, 
or Mb ~ —5.9 mag at 10 Gyr age (cf. Fig. Mb), corresponding to 
B = 25.9 at the distance of NGC 1316. As to disruption of clusters 
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Figure 15. _B-band luminosity function of the red cluster candidates in 
NGC 1316 after luminosity dimming to an age of 10 Gyr Hatching and 
hne types as in Fig. MSl 



Table 2. Evolution of luminosity and specific frequency of NGC 1316. 



Age 



2 Gyr 
10 Gyr 



AMv Sn 

Age(N1316) = 2Gyr* 



+0.00 
+1.63 



1.7 ±0.4 
8.2 ± 1.9 



AMv 

Age (Nl 316) = 



Sn 
2 Gyr + 14 Gyr* 



+0.00 
+0.14 



1.7 ± 0.4 
2.1 ±0.5 



Note to TaWeg: * The values listed under the column header 'Age (N1316) 
= 2 Gyr' were calculated under the assumption that all integrated light of 
NGC 1316 is 2 Gyr old, while the values listed under 'Age (N1316) = 2 Gyr 
+ 14 Gyr' were calculated under the assumption that only the integrated 
light within a galactocentric radius of 20" (which is currently 15% of the 
total light) is 2 Gyr old, while the rest is 14 Gyr old. 



during the evolution of a GC system, we refer to the simulations 
of Vesperini (2000) who showed that the great majority of the dis- 
rupted clusters are those with initial masses ^ 10^ Mq (cf. Fig. 13 
of Vesperini 2000), corresponding to Ms J; —5.2 mag for solar- 
metallicity GCs at an age of 10 Gyr, which is beyond the turnover. 
Hence, this effect should not alter the LF of the clusters that are 
brighter than the turnover significantly. 

Considering the above, we estimated the Sn of the GC system 
of NGC 1316 at a merger age of 10 Gyr as follows. For the 'blue' 
population, we counted the GCs up to the absolute magnitude of 
the turnover of the LF of our Galaxy (Mb = —6.6, i.e., B = 25.2), 
corrected for incompleteness as a function of magnitude, and mul- 
tiplied that number by a factor 2 (i.e., assuming a log-normal LF for 
the evolved GC system). The result is a total of 1950 ± 140 'blue' 
clusters. For the 'red' population, we counted the clusters having 
a 'faded' B magnitude up to 25.9 (as discussed above), and multi- 
plied that number by a factor 2 as well, leading to a total of 1450 
± 145 clusters. The corresponding Sjv values are listed in Table Q, 
which shows that even if the currently intermediate-age stellar pop- 
ulation in NGC 1316 only inhabits its central part, the specific clus- 
ter frequency will evolve to a value consistent with those of typical 
ellipticals in the field and poor groups (for which Sn = 2.6 ± 0.5, 
Harris 1991) at an age of ~ 10 Gyr. 



5 SUMMARY 

Using a combination of HST and large-field ground-based imaging 
of NGC 1316, a well-known early-type merger remnant galaxy for 
which we determined an age of -^ 3 Gyr from spectra of its brightest 
GCs in Paper 1, we presented the optical and near-IR photometric 
properties of its star cluster system in this paper. Our main results 
are itemized below, followed by a short discussion. 

• The B—I colour distribution of the clusters exhibits two peaks: 
One at 5 — / ~ 1 .5 (coinciding with the peak B~I colour of the 
Galactic GC system) and one at B — 7 ~ 1.8. The brightest clusters 
mainly populate the red peak, while the blue peak gets more popu- 
lated going towards fainter magnitudes. The clusters in the red peak 
make up about 50-60 per cent of the total cluster population. 

• By comparing the photometric properties of the brightest 8 
clusters with the those of the Galactic GC system as well as with 
current SSP models which include the contribution of the AGB 
phase of stellar evolution, we find that 

(i) They are up to an order of magnitude more luminous than 
ui Cen, the most massive GC in our Galaxy, and hence are most 
probably significantly younger than 10-14 Gyr; 

(ii) Their optical-near-lR colours are not as red as intermediate- 
age (0.3-1.5 Gyr old) GCs such as found in the Magellanic 
Clouds and in NGC 7252 (which are due to the presence of TP- 
AGB stars). Instead, their optical and near-lR colours are all 
best fit by a 2.5-3 Gyr old population of roughly solar metallic- 
ity, consistent with the results from our earlier spectral analysis 
(cf. Paper I). 

Together with the optical properties of the cluster system described 
in the previous item, these results are consistent with a scenario 
involving the presence of two cluster subpopulations in NGC 1316: 
(/) A 'red' (B~I ~ 1.8) population of clusters which were created 
during a major gas-rich merger that occurred ~ 3 Gyr ago, having 
solar metallicity to within ± 0.15-0.2 dex (as derived from spectra 
in Paper I), and (ii) a 'blue' (B~I ~ 1 .5) population of 'old' , metal- 
poor clusters that were associated with the halos of the pre-merger 
galaxies. We show that under this assumption, the 'red' clusters 
will evolve to B~I ~ 2.15 in a time span of 10 Gyr, which colour 
is remarkably similar to those of the red peaks observed in bimodal 
distributions of cluster colours in well-studied 'normal, old' giant 
ellipticals with galaxy luminosities similar to that of NGC 1316, 
suggesting that the red clusters in giant ellipticals have a similar 
origin. 

• The luminosity function (LF) of the cluster population is best 
represented by a power law with a slope a — — 1.2 ± 0.3. The 
power-law behaviour is dominated by the 'red' clusters; the LF of 
clusters populating the blue end of the colour distribution shows ev- 
idence for a 'turnover' as found in the GC systems of nearby galax- 
ies like the Milky Way and M 31. The power-law slope of the GC 
system ofNGC 1316 is flatter than those found in 'young' mergers 
with ages ^ 500 Myr. We suggest that this flattening of the power- 
law LF slope represents a stage in the dynamical evolution of star 
cluster mass functions that is intermediate between those found in 
'young' merger remnants and those found in 'normal', old giant el- 
lipticals (which show a log-normal LF). Using simple scaling rela- 
tions between time-scales of long-term cluster disruption processes 
(tidal shocking and evaporation), we estimate the absolute magni- 
tude of the turnover of the LF of a 3 Gyr old cluster population to 
be in the range —5.1 ^ Mb ^ —4.4, which is beyond the 50% 
completeness limit of the available WFPC2 data. 

• The radial distribution of the surface number density of 
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NGC 1316's star cluster system is consistent with that of its inte- 
grated light in the outer regions, while it flattens off slightly towards 
the centre relative to the integrated light. The latter flattening to- 
wards the centre is significantly less pronounced than in 'old' giant 
ellipticals, in which the outer surface density profile of clusters is 
typically more extended than that of the underlying starlight. These 
findings suggest that the cluster system of NGC 1316 (and, by anal- 
ogy, of giant ellipticals in general), originally experienced the same 
violent relaxation as did the main body of the merger remnant, af- 
ter which tidal shocking of (mosty inner) clusters cause the cluster 
surface density profile to become gradually flatter with time. 

• The present-day specific frequency of clusters with B < 26.0 
is Sn = 1.7 ± 0.4. After taking into account age dimming of the 
3-Gyr-old clusters as well as of the integrated light of the galaxy 
to an age of 10 Gyr [under the conservative assumption that only 
the inner 20" radius of the stellar population of NGC 1316 has 
a luminosity-weighted age of 2 Gyr (as derived spectroscopically 
by Kuntschner 2000) and that the rest of the light is from 14 Gyr 
old stars], the specific frequency will increase to 2.1 ± 0.5, which 
is consistent with those found in ellipticals in the field and poor 
groups. 

All these features of the star cluster system of NGC 1316 
constitute important evidence in support of 'merger scenarios' for 
forming giant elliptical galaxies through major gas-rich mergers. 
I.e., the system of luminous, metal-rich clusters in NGC 1316 is 
consistent with being a recently-formed analog of the 'red' clusters 
in 'normal' giant elliptical galaxies. 



is performed in situ by type II supemovae, which are the main pro- 
ducers of a elements. Hence, if major gas-rich mergers are indeed 
responsible for forming giant elliptical galaxies (and the second- 
generation clusters), those mergers most likely occurred up to a few 
Gyr after the pre-merger galaxies were created. This allowed time 
for the latter's ISM to become sufficiently enriched (most proba- 
bly primarily by SNII ejecta, cf. Thomas et al. 1999). New obser- 
vations should be conducted to find out whether or not the 'red' 
clusters in 'normal' giant ellipticals show supersolar a element 
abundances. If found, this will constitute an important constraint to 
merger models such as the well-known 'hierarchical' merger model 
for elliptical galaxies (e.g., Kauffmann 1996) in the sense that the 
building blocks that form giant elliptical galaxies must have already 
existed ~ 10^ yr after the 'halo' GCs around galaxies were formed. 
However, merging of two sizable gas-rich spirals at late 
epochs (which is relevant to young and intermediate-age merger 
remnants such as NGC 1 3 16 and NGC 7252) is nof expected to pro- 
duce supersolar Q-element abundances, because most of the stars in 
the merging spiral disks have formed in a long-lasting (~ 10 Gyr) 
star formation process, leading to approximately solar abundances 
and element ratios in their ISM. If stars form out of such gas with 
a Salpeter IMF, they are expected to show roughly solar [a/Fe] ra- 
tios. In such mergers at late epochs, Q-enhancement only occurs 
in models featuring a significantly flattened IMF during the major 
starburst following the merger (Thomas et al. 1999). At an age of 
~ 3 Gyr, the spectral features of Magnesium and Iron in the blue- 
visual region should be easily detectable. Thus, the bright clusters 
in NGC 1316 are very interesting probes to test such scenarios for 
chemical enrichment of merging galaxies. 



6 OUTSTANDING ISSUES 

The presence of very luminous second-generation GCs in 
NGC 1316 that are of an age at which dynamical evolution is likely 
to have already had a significant impact is relevant to a number of 
outstanding issues of interest going beyond the particular case of 
NGC 1316. One such issue is whether the mass function of clus- 
ters created in mergers is biased towards higher masses than that 
of clusters in 'normal' galaxies. E.g., W3, the brightest cluster in 
NGC 7252 (Whitmore et al. 1993) has a reported absolute magni- 
tude M^ = -16.2 (Miller et al. 1997). At i ts age of ~ 300 Myr and 
its metallicity of 0.5 Zq, Maraston et al. (2001) derive a mass of 
~ 3.7 X 10^ Mq (for a Salpeter IMF), which is ~ 10 times higher 
than the mass of uj Cen. Even for a significantly older merger rem- 
nant such as NGC I3I6, the brightest cluster (designated # 114 in 
Paper I) has My = —13.0, which translates to a mass of (1.4 ± 
0.2) X 10'^ Mq at an age of 3.0 ± 0.5 Gyr and solar metallic- 
ity (for a Salpeter IMF): still a factor ~ 5 higher than the mass of 
UJ Cen. The masses of these bright clusters will be checked by high- 
dispersion dynamical measurements using the ESO VLT in order to 
enable further study of mass functions of GCs created in mergers. 
This will be the subject of a forthcoming paper. 

Another important outstanding issue is whether or not such 
second-generation clusters show supersolar [a/Fe] element ratios. 
The stellar populations of many 'normal, old' giant ellipticals are 
known to be a-enhanced by ~ 0.2-0.3 dex (Worthey, Faber & 
Gonzalez 1992; Davies, Sadler & Peletier 1993; Kuntschner 2000). 
Recent models of galaxy formation that incorporate the evolution 
of chemical enrichment of the ISM during bursts (Thomas, Greg- 
gio & Bender 1999) show that such a-element enhancements can 
be reached by merging of small sub-galactic entities on a short time 
scale (^ 1 Gyr) so that much of the chemical enrichment of the ISM 



Acknowledgments. 

This work is based on observations obtained at the European South- 
em Observatory, La Silla, Chile (Observing Programme 58. E- 
0666), and on archival data of the NAS A/ESA Hubble Space Tele- 
scope, which is operated by AURA, Inc., under NASA contract 
NAS 5-26555. We are grateful to the European Southern Obser- 
vatory for allocating observing time to this project, and for having 
done a great j ob of archiving observations made with the NTT since 
its inauguration. We thank Chris Lidman and the technical support 
staff of the European Southern Observatory for their support dur- 
ing the near-IR observations. Brad Whitmore and Markus Kissler- 
Patig for very useful discussions, and the referee Jean Brodie for 
comments that led to an improvement of the paper. MVA was 
partially supported by grants from CONICET, Agenda Nacional 
de Promocion Cientifica y Tecnologica and Secretaria de Cien- 
cia y Tecnica de la Universidad Nacional de Cordoba, Argentina. 
MVA also thanks the Chilean Fondecyt for a visitor grant. CM 
is supported by the 'Sonderforschungsbereich 375-95 fiir Astro- 
Teilchenphysik' of the Deutsche Forschungsgemeinschaft. DM's 
work was performed in part under the auspices of the Chilean 
Fondecyt Nos. 01990440 and 07990048, and of DIPUC No. 98.16e. 



REFERENCES 

Aaronson, M., Cohen, J. G., Mould, J., Malkan, M., 1978, ApJ, 223, 824 
Aguilar, L. A., White, S. D. M., 1985, ApJ, 295, 374 
Aguilar, L. A., Hut, P., Ostriker, J. P., 1988, ApJ, 335, 720 
Amaboldi, M., Freeman, K. C, Gerhard, O., et al., 1998, ApJ, 507, 759 
Ashman, K. M., Zepf, S. E., 1992, ApJ, 384, 50 

Ashman, K. M., Zepf, S. E., 1998, Globular Cluster Systems. Cambridge 
University Press, Cambridge 



© 0000 RAS, MNRAS 000, 000-000 



Star clusters in NGC 1316 19 



Brodie, J. P., Schroder, L. L., Huchra, J. P., Phillips, A. C, Kissler-Patig, 

M., Forbes, D. A., 1998, AJ, 116, 691 
Bruzual, G. A., Chariot, S., 1993, ApJ, 405, 538 
Bruzual, G. A., Chariot, S., 1996, private communication (BC96) 
Burstein, D., Heiles, C, 1984, ApJS, 54, 33 
Caon, N., Capaccioh, M., D'Onofrio, M., 1994, A&AS, 106, 199 
Carter, B. S., Meadows, V. S., 1995, MNRAS, 276, 734 
Cohen, J. G., 1982, ApJ, 258, 143 

Cohen, J. G., Frogel, J. A., Persson, S. E., Elias, J. H., 1981, ApJ, 249, 481 
Cote, P, Marzke, R. O., West, M. J., 1998, ApJ, 501, 554 
Davies, R. L., Sadler, E. M., Peletier, R. R, 1993, MNRAS, 262, 650 
D'Onofrio, M., Capaccioli, M., Zaggia, S. R., Caon, N., 1997, MNRAS, 

289, 847 
Drinkwater, M. J., Gregg, M. D., Colless, M., 2001, ApJ, 548, L139 
Elmegreen, B. G., Efremov, Yu. N., 1997, ApJ, 480, 235 
Elson, R. A. W., Fall, S. M., 1985, PASP 97, 692 
Elson, R. A. W., Santiago, B., 1996, MNRAS, 278, 617 
Elson, R. A. W., Grillmair, C. J., Forbes, D. A., Rabban, M., Williger, G. 

M., Brodie, J. R, 1998, MNRAS, 295, 240 
Fall, S. M., Rees, M. J., 1977, MNRAS, 181, 37 
Ferguson, H. C, 1989, AJ, 98, 367 

Ferrarese, L., Ford, H. C, Huchra, J. P, et al., 2000, ApJS, 128, 431 
Fleming, D. E. B., Hanis, W. E., Pritchet, C. J., Hanes, D. A., 1995, AJ, 

109, 1044 
Forbes, D. A., Forte, J. C, 2001, MNRAS, 322, 257 
Forbes, D. A., Brodie, J. P, Grillmair, C. J., 1997, AJ, 113, 1652 
Forbes, D. A., Grillmair, C. J., Williger, G. M., Elson, R. A. W., Brodie, J. 

P, 1998, MNRAS, 293, 325 
Forte, J. C. Martinez, R. E., Muzzio, J. C, 1982, AJ, 87, 1465 
Frogel, J. A., Mould, J. R., Blanco, V. M., 1990, ApJ, 352, 96 
Gallagher, J. S. Ill, Smith, L. J., 1999, MNRAS, 304, 540 
Geisler, D., Lee, M. G., Kim, E., 1996, AJ, 111, 1529 
Glass, I. S., 1984, MNRAS, 211, 461 
Gnedin, O. Y., Ostriker, J. R, 1997, ApJ, 474, 223 
Goudfrooij, P., Hansen, L., J0rgensen, H. E., N0rgaard-nielsen, H. U., de 

Jong. T., van den Hoek, L. B., 1994a, A&AS, 104, 179 
Goudfrooij, P., Mack, J., Kissler-Patig, M., Meylan, G., Minnid, D., 2001, 

MNRAS, 643, 657 (Paper I) 
Grillmair, C. J., Faber, S. M., Lauer, T. R., et al., 1994, AJ, 108, 102 
Grillmair, C. J., Forbes, D. A., Brodie, J. P, Elson, R. A. W., 1999, AJ, 1 17, 

167 
Harris, W. E.. 1991, ARA&A, 29, 543 
Harris, W. E., 1996, AJ, 112, 1487 
Harris, W. E., van den Bergh, S., 1981, AJ 86, 1627 
Hibbard, J. E., Mihos, J. C, 1995, AJ, 110, 140 
Holtzman, J. A.. Faber, S. M., Shaya, E. J., et al. 1992, AJ, 103, 691 
Holtzman, J. A., Burrows, C. J., Casertano, S., et al., 1995, PASP 107, 165 
Iben, I. Jr, Renzini, A., 1983, ARA&A, 21, 271 
Jog, C. J., Solomon, P M., 1992, ApJ, 387, 152 
Kauffmann, G., 1996, MNRAS, 281, 487 
Kundu, A., Whitmore, B. C, Sparks, W. B., et al 1999. ApJ. 513 733 



Kundu, A., Whitmore, B. C, 2001, AJ, in press (^stro-ph/0103021|) 

Kuntschner, H., 2000, MNRAS, 315, 184 

Langon, A., 1999, in Le Bertre T., Lebre A., and Waelkens C, Asymptotic 
Giant Branch Stars. Kluwer, Dordrecht, p. 579 

Landolt, A. U., 1992, AJ, 104, 340 

Larsen, S. S., Brodie, J. P.. Huchra. J. P .. Forbes, D. A., Grillmair, C. J., 
2001, AJ, in press ( sistro-ph/0 102374 ) 

Maraston, C, 1998, MNRAS, 300, 872 

Maraston, C, 1999, in Carral, P. & J. Cepa, eds, Star formation in Early- 
type Galaxies. ASP, San Francisco, p. 28 

Maraston, C, Kissler-Padg, M., Brodie, J. P., et al., 2001, A&A, in press 

Meurer, G. R., Heckman, T. M., Leitherer, C, et al., 1995, AJ, 110, 2665 

Meylan, G., Heggie, D. C, 1997, A&AR, 8, 1 

Miller, B. W., Whitmore, B. C, Schweizer, F, Fall, S. M., 1997, AJ, 114, 
2381 

Minnid, D., 1995, AJ, 109, 1663 



Minniti, D., Alonso, M. V., Goudfrooij, P., Jablonka, P., Meylan, G., 1996, 

ApJ, 467, 221 
Murah, C, Weinberg, M. D., 1997, MNRAS, 288, 749 
Olszewski, E. W., 1995, in Gilmore G., van der Kruit P., eds., Stellar Popu- 
lations. Kluwer, Dordrecht, p. 181 
Persson, S. E., Aaronson, M., Cohen, J. G., Frogel, J. A., Mathews, K., 

1983, ApJ, 266, 105 
Renzini, A., 1981, Ann. Phys. Fr, 6, 87 
Renzini, A., Voli, M., 1981, A&A, 94, 175 
Renzini, A., Buzzoni, A., 1986, in Chiosi C, Renzini A., Spectral Evolution 

of Galaxies. Reidel, Dordrecht, p. 195 
Riess, A. G., Filippenko, A. V., Challis, P, et al., 1998, AJ, 1 16, 1009 
Salpeter, E. E., 1955, ApJ, 121, 161 
Sandage, A. R., Tammann, G. A., 1987, A Revised Shapley-Ames Catalog 

of Bright Galaxies. Carnegie Institution of Washington, Washington, 

D.C. 
Sanders, D. B., Mirabel, I. F, 1996, ARA&A, 34, 749 
Scalo, J. M., 1986, Fundam. Cosmic Phys., 11, 1 
Schweizer, F, 1980, ApJ, 237, 303 
Schweizer, R, 1981, ApJ, 246, 722 
Schweizer, F, 1987, in Faber S. M., Nearly Normal Galaxies. Springer, New 

York, p. 18 
Schweizer, F., 1998, in Friedli D., Martinet L., Pfenniger D., eds, 'Galaxies: 

Interactions and Induced Star Formation'. Springer, Berlin, p. 105 
Schweizer, F, Seitzer, P, 1998, AJ, 116, 2209 
Schweizer, F, Miller, B. W., Whitmore, B. C, Fall, S. M., 1996, AJ, 112, 

1839 
Searle, L., Wilkinson, A., Bagnuolo, W. G., 1980, ApJ, 239, 803 (SWB) 

Schlegel, D. J., Finkbeiner, D. P, Davis, M., 1998, ApJ, 500. 525 

Smith, L. J., Gallagher, J. S. Ill, 2001, MNRAS, in press ( ^stro-ph/0 104429 1 

Stetson, P B., 1987, PASR 99, 191 

Spitzer, L., Chevalier, R. A., 1973, ApJ, 183, 565 

Terievich, A. I, Kuntschner, H., Bower, R. G., Caldwell, N., Sharpies, R. 

M., 1999, MNRAS, 310, 445 
Thomas, D., Greggio, L., Bender, R., 1999, MNRAS, 302, 537 
van den Bergh, S., 1981, A&AS, 46, 79 
van den Bergh, S., 1995, ApJ, 450, 27 
Vesperini, E., 2000, MNRAS, 318, 841 
Wegner, G., 1979, Ap&SS, 60, 15 
Westerlund, B. E., Azzopardi, M., Rebeirot, E., Breysacher, J., 1991, 

A&AS, 91, 425 
Whitmore, B. C, 1997, in Livio M., ed.. The Extragalactic Distance Scale. 

Cambridge University Press, Cambridge, p. 254 
Whitmore, B. C, Schweizer, F, Leitherer, C, Borne, K., Robert, C, 1993, 

AJ, 106, 1354 
Whitmore, B. C, Miller, B. W., Schweizer, F, Fall, S. M., 1997, AJ, 114, 

1797 
Whitmore, B. C, Heyer, I., Casertano, S., 1999a, PASP 111, 1559 
Whitmore, B. C, Zhang, Q., Leitherer, C, Fall, S. M., Schweizer, F, Miller, 

B. W., 1999b, AJ, 118, 1551 
Worthey, G., Faber, S. M., Gonzalez, J. J., 1992, ApJ, 398, 69 
Zepf, S. E., Carter, D., Sharpies, R. M., Ashman, K. M., 1995a, ApJ, 445, 

L19 
Zepf, S. E., Ashman, K. M., Geisler, D., 1995b, ApJ, 443, 570 
Zepf, S. E., Ashman, K. M., English, J., Freeman, K. C, Sharpies, R. M., 

1999, AJ, 118,752 



APPENDIX A: PHOTOMETRY TABLES OF STAR 
CLUSTER CANDIDATES 



© 0000 RAS, MNRAS 000, 000-000 



20 Paul Goudfrooij et al. 

Table Al. Photometry and astrometry of the 50 brightest GC candidates on the WFPC2 frames. JHK photometry is from the 2.2-m ESO/MPI telescope. The 
object list is sorted on B magnitude (brightest first). 



ARA 


A DEC 


B 


B-I 


K 


J~K 


H-K 


arcsec 


arcsec 


mag 


mag 


mag 


mag 


mag 


8.91 


-11.46 


19.627 ± 0.003 


1.866 ±0.005 


15.78 ±0.04 


0.83 ± 0.06 


0.15 ±0.06 


-48.09 


-22.24 


20.241 ± 0.003 


1.744 ±0.004 


16.60 ± 0.06 


0.78 ± 0.08 


0.11 ±0.08 


-25.01 


-13.04 


20.343 ± 0.003 


1.700 ± 0.005 


16.75 ± 0.07 


0.95 ± 0.09 


0.18 ±0.09 


n.n 


12.14 


20.467 ± 0.003 


1.935 ±0.005 


15.83 ±0.05 


0.95 ± 0.08 


0.35 ± 0.07 


-43.94 


-38.32 


20.523 ± 0.003 


1.740 ±0.005 


16.73 ±0.06 


0.83 ± 0.08 


0.28 ± 0.08 


4.13 


10.50 


20.918 ± 0.005 


2.103 ±0.005 


16.74 ±0.07 


1.06 ±0.11 


0.22 ±0.10 


-45.15 


-21.03 


20.967 ± 0.005 


1.757 ±0.006 


17.41 ±0.10 


0.86 ±0.13 


0.23 ±0.13 


53.03 


23.58 


21.313 ±0.007 


1.972 ±0.008 








-12.77 


-3.75 


21.609 ± 0.008 


2.047 ± 0.010 








-59.81 


-33.74 


22.094 ± 0.008 


1.721 ±0.011 








3.99 


74.26 


22.373 ± 0.009 


1.959 ±0.011 








-15.88 


-9.49 


22.413 ±0.012 


1.618 ±0.017 








-49.09 


-31.98 


22.430 ± 0.010 


1.754 ±0.013 








1.02 


-13.66 


22.483 ± 0.009 


1.918 ±0.013 








37.20 


72.27 


22.654 ± 0.010 


1.740 ±0.013 








1.74 


33.59 


22.656 ±0.013 


1.909 ±0.015 








-63.76 


21.37 


22.657 ±0.011 


1.771 ±0.014 








20.40 


63.09 


22.706 ± 0.012 


2.178 ±0.013 








-25.66 


-9.95 


22.729 ± 0.016 


2.071 ±0.019 








2.52 


-9.22 


22.743 ± 0.015 


1.883 ±0.022 








43.14 


72.73 


22.770 ±0.011 


2.125 ±0.013 








-23.15 


-39.81 


22.870 ± 0.015 


1.813 ±0.019 








13.54 


18.60 


22.908 ± 0.023 


1.605 ±0.028 








-6.76 


-15.93 


22.950 ± 0.014 


2.430 ± 0.016 








-38.30 


64.33 


22.951 ±0.013 


1.469 ±0.017 








29.16 


46.12 


23.006 ± 0.016 


1.868 ±0.018 








-5.48 


-7.77 


23.030 ± 0.021 


1.916 ±0.025 








-41.55 


-23.01 


23.041 ± 0.017 


1.677 ±0.021 








-62.41 


-54.32 


23.066 ± 0.015 


1.604 ±0.018 








-24.13 


2.52 


23.074 ± 0.022 


1.781 ±0.026 








-24.49 


78.68 


23.089 ± 0.013 


1.534 ±0.017 








-8.07 


-7.71 


23.131 ±0.022 


2.207 ± 0.025 








-28.74 


-13.85 


23.134 ±0.022 


1.867 ±0.026 








-21.48 


-42.39 


23.197 ±0.021 


1.381 ±0.027 








-27.23 


-2.73 


23.256 ± 0.026 


1.493 ±0.035 








-52.60 


-13.63 


23.276 ± 0.020 


1.636 ±0.024 








15.78 


-6.93 


23.278 ± 0.023 


1.615 ±0.030 








-34.60 


-27.10 


23.288 ± 0.023 


1.754 ±0.027 








-21.84 


75.37 


23.292 ± 0.017 


1.644 ±0.021 








-0.27 


13.01 


23.295 ± 0.022 


2.309 ± 0.025 








-36.77 


5.88 


23.302 ± 0.021 


1.5 14 ±0.028 








-82.63 


-46.10 


23.302 ± 0.016 


1.642 ±0.020 








13.36 


19.33 


23.311 ±0.033 


1.958 ±0.039 








53.16 


79.19 


23.320 ±0.018 


1.328 ±0.023 








-41.75 


38.95 


23.327 ± 0.018 


1.674 ±0.023 








-85.71 


52.72 


23.350 ±0.015 


1.722 ±0.019 








-21.31 


13.12 


23.352 ± 0.026 


1.464 ±0.034 








12.11 


45.72 


23.379 ± 0.019 


1.230 ±0.028 








-54.35 


-52.15 


23.389 ± 0.023 


1.790 ±0.027 








-49.59 


7.15 


23.394 ± 0.019 


1.641 ±0.026 









Note to TaWeUfl: The first two columns give the positional offsets of the GC candidates from the centre of NGC 1316 in RA 
and DEC, respectively. 



© 0000 RAS, MNRAS 000, 000-000 



Star clusters in NGC 1316 21 



Table A2. Photometry and astrometry of the 50 brightest GC candidates on the ESO NTT images. GC candidates that were also present on the WFPC2 images 
have been excluded from this list. JHK photometry is from the 2.2-m ESO/MPl telescope. The object list is sorted on B magnitude (brightest first). 



ARA 


A DEC 


B 


B-V 


B-I K 


J-K H~K 


arcsec 


arcsec 


mag 


mag 


mag mag 


mag mag 


-158.09 


64.95 


20.17 ±0.19 


0.98 ± 0.20 


1.78 ±0.20 




50.75 


-46.88 


20.64 ± 0.23 


0.87 ± 0.25 


1.45 ±0.25 




-84.07 


-163.18 


20.71 ±0.24 


0.97 ± 0.26 


1.73 ±0.25 




148.33 


132.35 


20.79 ± 0.25 


0.90 ± 0.27 


1.67 ±0.26 




53.89 


23.29 


20.83 ± 0.01 


1.04 ±0.09 


1.82 ±0.01 




53.89 


23.29 


20.83 ± 0.26 


1.04 ±0.27 


1.82 ±0.27 




-60.97 


94.40 


20.88 ± 0.26 


0.93 ± 0.28 


1.70 ±0.28 




67.31 


54.96 


21.03 ±0.28 


0.78 ± 0.30 


1.37 ±0.30 




67.31 


54.96 


21.03 ±0.28 


0.78 ± 0.30 


1.37 ±0.30 




178.46 


52.70 


21.15 ±0.30 


0.72 ± 0.32 


1.13 ±0.32 




86.04 


191.25 


21.22 ±0.31 


0.72 ± 0.33 


1.01 ±0.34 




160.04 


-150.52 


21.23 ±0.31 


0.64 ± 0.34 


0.99 ± 0.34 




-123.07 


8.43 


21.50 ±0.35 


1.01 ±0.37 


1.83 ±0.36 




40.89 


-198.75 


21.57 ±0.36 


0.68 ± 0.39 


1.02 ±0.39 




-141.71 


-0.40 


21.60 ±0.36 


1.20 ±0.38 


2.06 ± 0.38 




12.87 


129.47 


21.71 ±0.38 


0.86 ± 0.41 


1.51 ±0.41 




203.34 


21.51 


21.73 ±0.39 


0.72 ± 0.42 


0.99 ± 0.43 




71.17 


20.14 


21.75 ±0.39 


0.85 ± 0.42 


1.64 ±0.41 




-127.28 


-31.86 


21.79 ±0.40 


0.97 ± 0.42 


1.71 ±0.42 




60.05 


17.62 


21.83 ±0.41 


0.97 ± 0.44 


1.71 ±0.43 




5.25 


74.90 


21.96 ±0.01 


0.85 ±0.17 


1.57 ±0.01 




129.89 


-142.12 


21.97 ±0.44 


0.31 ±0.49 


1.23 ±0.47 




93.83 


49.74 


21.97 ±0.43 


0.92 ± 0.46 


1.71 ±0.46 




-171.13 


-176.59 


22.00 ± 0.44 


1.09 ±0.46 


1.99 ±0.46 




-15.98 


-49.92 


22.07 ± 0.45 


0.99 ± 0.48 


1.86 ±0.47 




-140.28 


58.05 


22.08 ± 0.46 


0.99 ± 0.49 


1.73 ±0.48 




-9.71 


-49.94 


22.10 ±0.46 


0.96 ± 0.49 


1.82 ±0.48 18.08 ±0.15 


1.23 ± 0.20 0.50 ± 0.20 


94.24 


5.95 


22.15 ±0.47 


1.14 ±0.50 


2.02 ± 0.49 




-133.37 


-106.62 


22.16 ±0.47 


0.87 ±0.51 


1.52 ±0.50 




-113.18 


80.14 


22.17 ±0.47 


0.90 ±0.51 


1.62 ±0.50 




190.76 


9.15 


22.21 ± 0.49 


0.56 ± 0.53 


1.38 ±0.52 




-119.69 


-7.11 


22.34 ±0.51 


1.02 ±0.55 


1.78 ±0.54 




-76.75 


-85.49 


22.50 ± 0.55 


1.02 ±0.59 


1.59 ±0.58 




69.04 


-6.28 


22.51 ±0.56 


0.84 ± 0.60 


1.58 ±0.59 




-115.13 


112.73 


22.51 ±0.56 


1.01 ±0.59 


1.79 ±0.58 




-120.31 


-70.95 


22.55 ± 0.57 


0.91 ±0.61 


1.69 ±0.60 




-93.24 


-2.41 


22.57 ± 0.57 


1.01 ±0.61 


1.77 ±0.60 




-100.24 


77.63 


22.59 ± 0.57 


1.05 ±0.61 


1.91 ±0.60 




-110.05 


-148.12 


22.60 ± 0.58 


1.07 ±0.61 


1.82 ±0.61 




-156.04 


-17.66 


22.60 ± 0.58 


0.79 ± 0.63 


1.34 ±0.62 




122.96 


82.64 


22.61 ± 0.59 


0.61 ±0.64 


1.22 ±0.63 




-33.49 


-164.70 


22.63 ± 0.59 


0.91 ± 0.63 


1.64 ±0.62 




-37.41 


-59.81 


22.63 ± 0.59 


0.99 ± 0.63 


1.54 ±0.63 




110.55 


1.67 


22.67 ± 0.60 


1.08 ± 0.63 


1.72 ±0.63 




-18.57 


108.55 


22.67 ± 0.60 


0.86 ± 0.64 


1.72 ±0.63 




-10.99 


-63.23 


22.73 ± 0.63 


0.63 ± 0.68 


1.15 ±0.68 




-104.57 


-106.84 


22.75 ± 0.62 


0.41 ± 0.69 


1.03 ±0.68 




100.28 


56.20 


22.75 ± 0.63 


0.57 ± 0.68 


1.32 ±0.67 




140.89 


114.79 


22.77 ± 0.63 


0.74 ± 0.68 


1.39 ±0.67 




27.38 


-104.95 


22.81 ±0.64 


0.96 ± 0.68 


1.68 ±0.67 





Note to Table |A2| The first two columns give the positional offsets of the GC candidates from the centre of NGC 1316 in RA and DEC, 
respectively. 
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